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1. INTRODUCTION 


The purpose of this program, BIOMAP, is to facilitate the 
electron microprobe analyses of biological specimens by 
computing from microprobe intensity count data and other 
parameters the quantity of specified chemical elements con- 
tained in formed elements of the blood and in micro-thin 
sections of tissue and gelatin. In addition, BIOMAP esti- 
mates background and concentration intensity counts for 
chemical elements in standards and computes and plots the 
distributions of elemental concentrations in red blood cell 
s amp 1 e s . 

BIOMAP was developed as part of the Preventative Medicine 
Division Project 1030, Intracellular Chemical Characteriza- 
tion, which calls for research and development in areas of 
cellular chemical analysis by micro-X-ray spectrometry and 
UV mi crospectrometry . Direction, of this effort was provided- 
by Dr. Stephen L. Kimzey, Chief of the Cellular Analytical 
Section of the Clinical Laboratory, Preventative Medicine 
Division, NASA-MSC. Electron microprobe data used to test 
the mathematical model were provided by Mrs. Linda Burns, 
electron microscopist , The University of Texas Medical School 
at Galveston. The physical analysis, mathematical modeling, 
and programming were done by the author. 

The atomic constant data and the Fortran-coded algorithm 
for calculating mass absorption coefficients are from the 
quantitative electron microprobe analysis program, MAGIC, 
written by J. W. Colby of Bell Telephone Laboratories , Inc., 
Allentoxirn, Pennsylvania. 


$AGE BLANK NOT 
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2. PROGRAM DESCRIPTION 


2. 1 General Description 

The central feature of BIOMAP is the computation of 
absorption factors and the iterative computation of unknown 
elemental concentrations. Absorption factors quantify the 
amounts of characteristic X-radiation absorbed by the sample 
or standard and correctly relate X-ray intensity counts with 
elemental concentrations in the microprobe equation. Com- 
putation of unknown elemental concentrations is accomplished 
by solving the quantitative equation. This equation is 
functionally similar to the microprobe equation but involves 
instead ratios of the intensity counts, absorption factors, 
and elemental concentrations from unknown and standard 
samples . 

Since a significant degree of variation is characteristic 
of both input and output data, certain statistical features 
have also been built into the program. First is a general 
linear hypothesis subroutine which inputs intensity count 
and absorption corrected area density data and computes 
estimates of the residual background and the rate of count 
increase with area density, along with associated variances. 
From these estimates and variances an estimate of the most 
probable intensity count for a given standard and its asso- 
ciated variance can be made. Second, if the unknown sample 
analyzed is a set of red blood cells (RBC) j certain statis- 
tical parameters of the sample elemental concentrations are 
computed and cumulative and frequency functions of these 
concentrations are plotted. 
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BIOMAP can be set up to provide the following types of 

individual or combination analyses: 

(1) Analyze microprobe data from a set of standards to 
statistically estimate the residual background intensity 
count and the concentration intensity count for the 
analyzed elements. 

(2) Analyze microprobe data from a set of gelatin samples 
of known or unknown chemical content, obtaining con- 
centration of specified elements. 

(3) Analyze microprobe data from a set of red blood cells, 
obtaining weight percentages and weight percentage 
distribution plots. 

Analyses (2) and (3) depend on input data from a previous 

type (1) analysis; hoxvever, combination analyses of types 

(1) and (2) and of types (1) and (3) can be made. 

The logical sequence of program operations is as follows: 

(1) References atomic constant data from block data storage. 

(2) Provides by data statement the weight fractions of 
chemical elements in normal red blood cells and plain 
gelatin 

(3) Reads data cards to input problem information, micro- 
probe parameters, sample and standard data, and program 
mode 

(4) Constructs atomic number arrays for analyzed elements 
and elements added to standards 

(5) Forms symbol, atomic number, and area density arrays 
for elements in dry gelatin standards 
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(6) Computes the mass absorption coefficients of the 
sample and standard elements for the X-ray wavelengths 
analyzed 

(7) Computes partial and total mass absorption coefficients 
of standard for each analyzed X-ray wavelength 

(8) Computes absorption factors for each analyzed X-ray 
wavelength of standard 

(9) Computes area densities of analyzed elements in 
standard 

(10) Computes estimates of statistical parameters for a 
set of standards and the best linear unbiased esti- 
mate of the accumulated count for a specified standard. 

For analysis of red blood cells and thin gelatin or tissue 

samples , the program performs the following operations : 

(11) Inputs chemical symbols and intensity counts for 
analyzed elements in each sample 

(12) Forms initial area density arrays for each sample 

(13) Computes mass absorption coefficient of sample for 
lines analyzed 

(14) Computes absorption factors for analyzed lines in sample 

(15) Computes the absorption correction for each element 
analyzed 

(16) Computes the area density of each analyzed element in 
the sample 

(17) Updates area density and tests against reiteration 
criterion 

(18) Computes the volume density of analyzed elements in 
gelatin samples 
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(IS) Computes weight fractions of analyzed elements in red 
blood cell sample set 

(20) Computes weight -percentage range, average, variance, 
and percentage standard deviation over red blood cell 
set for each analyzed element 

(21) Provides cumulative distribution and histogram plots 
of the weight-percentage distributions. 

In operation (1) , the atomic constants read from thO'data 
tape for each element are the atomic number, the elemental 
symbol, the critical excitation potentials for the K, L, 
and M X-ray states, the wavelengths of the Ka, La, and Ma 
spectral lines, the absorption edge wavelengths, exponents, 
and empirical constants for mass absorption coefficient 
equations, and fifth-degree polynomial coefficients for the 
backscatter factor. 

In operation (3) , problem information consists of alpha- 
numerics for problem type, problem number, date, submitter's 
name, and a descriptive title. Microprobe parameters are 
accelerating potential', beam diameter, beam current, count- 
ing time, and counter dead time. Sample information con- 
sists of a sample type index, the number of .amples to be 
analyzed, the number of analyzed elements in each sample, 
and the wet thickness of the samples in microns. Standard 
information consists of the grams of gelatin per liter used 
in preparing the standards , the water fraction in the 
gelatin, the thickness of the wet standard sections, the 
number of standards input for statistical evaluation, and 
the symbols of the analyzed elements. If a single standard 
is used the estimated residual background and intensity 
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count for each standard element is read in with the standard 
symbols. The program parameters specify whether to print out 
intermediate results in namelists and the weight options in 
subroutine GLHFR2. If multiple standards are analyzed to 
obtain a linear unbiased estimate of the residual background 
and intensity count for each standard element, the number of 
elements added to each standard, and their chemical symbol ana 
amounts in grams are read in at the top of the standard 
analysis loop. Also, at the top of this loop is read the flag 
for the primary standard and the chemical symbols of the ana- 
lyzed standard elements and their concentration intensity 
count. If multiple standards are processed to estimate the 
primary standard intensity count, microfilm plots are produced 
of the data points and the fitted line. 

2 . 2 Technical Description 

2.2.1 Analysis 

In general, X-ray states are excited in microprobe specimens 
primarily by incident electrons and secondarily by suffi- 
ciently energetic X-ray photons resulting from the fluorescent 
decay of X-ray states. The average number n of specified 
X-ray states excited in a given atomic species per single- 
incident electron can be represented as 

n = n* + n M (1) 

where n T is the average number resulting from primary 
excitation and n" is the average number resulting from 
secondary excitation. 


2-5 



Equation (1) can be expressed in the alternate form 


n = n'f (2) 

where f > 0 is defined as the fluorescence factor. 

The derivation of an expression for n' is of central 
importance in quantitative microprobe theory and proceeds 
as follows: 

In differential form, the average number of specified X-ray 

states dn 1 of energy E excited within an electron path 

c 

length dx is 


dn ’ = <$Q(E(x),E c ) dx (3) 

where E is the average electron energy over dx; Q is 
the ionization cross-section of the given atomic X-ray 
state, a function of E and E ; and 6 is the number 

G 

of specified atoms per unit volume of specimen. 

In a thick specimen where the energy of the penetrating 
electron may decay below 3E c tfhile traversing a path 
length x , Q cannot be considered constant, and 

Q dx = 6 j [Q/CdE/dx)] dE 

s o 

For a thin low-density specimen, however, one can reasonably 
assume that entering electrons completely penetrate the 
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target with little loss of energy or scattering. As a 
consequence, Q is effectively constant and the path 
length dx can be equated to the sample thickness z . 

These approximations enable writing Equation (3) as 

n T = SQz (4) 

The average number of states produced in time t by a 
beam current i , considering a fraction of ionizations 
(1 - R) to be lost as a result of electron backscattering, 
can be expressed as 


N 1 = 4>tR6Qz (5) 

where <j> the electron flux equals i divided by the 
electronic charge q 

0 

Since the microprobe measures ionizing events caused by 
X-ray photons resulting from the decay of the excited 
X-ray states, the microprobe equation relating measured 
intensity and atomic concentration is 

I = (f2/4iT)Dwpc!>tafR<5Qz (6) 

These additional factors are defined as follows : 

(£2/4tt) is the fraction of the total solid angle intercepted 
by the X-ray spectrometer, D is the spectrometer-detector 
efficiency, w is the fluorescence yield of the particular 
X-ray state, and p is the probability that the photon 
emitted will be that specified. The factors a and f are 
the absorption and fluorescence factors respectively and both 
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depend on the chemical and physical configuration of the 
specimen. The absorption factor a is the probability of 
a photon not being absorbed as it passes through the sample, 
and is less than one. The fluorescence factor f accounts 
for the number of states produced secondarily by X-ray 
photons with energy greater than E , and is greater than 
one . 

A thin standard sample having a known concentration of the 
specified chemical element would yield a measured intensity 

I 0 - C^/4v)Di^ppta 0 f 0 R 0 z 0 Q6 0 (7) 

By forming the ratios I/I the common factors may be elimi- 
nated by division, leaving 

I Q /I = (a 0 /a) Cf 0 /fHV R} ( z o /z ) V 6 W 

If the standard is chemically constituted sufficiently 
similar to the unknown sample, the value of the factors 
(f /f) and (R q /R) , known as the fluorescence and backscatter 
corrections respectively, become essentially one, leaving 

V 1 - O 0 /a)(z 0 /z) <S 0 /6 (9) 

Gelatin-based standards can be chemically constituted 
similar to a biological specimen but cannot be cut as thin 
as a red blood cell for instance and, therefore, expres- 
sions must be derived for the absorption factor which 
depends on sample thickness. 
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An expression for the absorption factor can be simply derived 
by assuming that X-ray states are produced uniformly over the 
thickness z , a reasonable assumption since Q is essen- 
tially constant for E > 3E c « One can start with the differ- 
ential relation 

dl = qN • dz (10) 

when dl ' are the number of photons observed by the detector, 

N is the number of photons generated per unit path length, 
and q is the probability of a photon not being absorbed 
or scattered. 

The probability q can be identified from Beer’s law 

-XPs 

x = i Q e 

— y n 2 

as the exponential factor e where 

X = esc w i u i (11) 

and is the angle of emergence of the photons, is 

the weight fraction of the i th element in the sample and 
p. is the mass absorption coefficient of the i th element 
for the specified X-ray photon. 

By "integrating over the thickness z the number of 
photons detected I are 

I = N J e" Xps ’ dz 1 = N (1 - e“ Xpz )/xP 

o 


r 
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and the fraction detected (not absorbed) is 

a * Wz * ^ “ e" xpz )/xpz 
The absorption correction can be expressed as 


Ca Q /a) 


1 - e 


~X A P„z 


oo) 


Xpz 


(i 


*Xpz 


) x o p o z o 


( 12 ) 


which, when substituted into the quantitative equation, 
Equation (9) , yields 



(l - e~ W °)( XP z) Vo 
(1 * e- XPZ )(Wo) ' 6Z 


(13) 


Since biological specimens are prepared wet and dry under 
vacuum, it is no longer possible to know the resulting 
thickness or dry-volume density of a standard element. 
However, if there is no lateral redistribution of the sample 
material, the dry area densities = pgZ Q of the standard 

elements remain invariant and Equation (13) can be written 
and solved in terms of area densities 


(l - e X ° °)x :L q a i 
(> - ° 


(14) 


Here a and c 0 without superscripts represent the total 
area density of sample and standard respectively, and the 
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superscript i designates quantities specific for the i th 
element. 

2.2.2 Method of Solution 

If one knows as a zero order approximation the normal 
average-weight fractions of elements in a biological sample, 
Its total area density, and can analyze simultaneously any 
elements whose concentrations deviate significantly from 
normal, then the quantitative microprobe equation used with 
an iterative computational procedure will provide a higher 
o?der estimate of the area densities of the analyzed ele- 
ments. The iterative procedure is made necessary by the 
nature of the absorption coefficient which cannot be com- 
puted exactly before the solution is obtained. Convergence 
to solutions lying within the precision of experimental 
measurement is rapidly achieved, however, by recomputing the 
absorption factor of the sample subsequent to each iterative 
area density computation, 

A summary of the input data and computational method in 
mathematical notation is as follows: 

I. Computations related to the standard 
A. Input data 

1. w. —weight fractions of elements in gelatin 

(i = 1, L) 

2. h — fraction of v?ater in gelatin 

3. g — grams gelatin per liter standard 

solution 
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4. t 


— thickness of wet gelatin slice 


5. a. ~ gram atoms of standard elements added 

(j - 1, .... M) 

6. B k — accumulated background intensity count for 

k th analyzed element (k = 1, . , . , N) 

7 . N k — accumulated intensity count for k th 

analyzed element (k = 1 , . . . , N] 

8. Atomic constants ( J . W. Colby] 

B. Computations 

1. Compute mass of each element in gelatin, 

m!^ = vj\(1 - h)g (i = 1,...,L) 

2. Compute mass of each added element, 

mu = a. A. (j “ 1 , . . . ,M) 

3 3 3 

3. Compute total mass of each element in liter of 
standard solution. 

m. = m! + m'. 1 

X XX 

4. Compute volume density of each element in 
standard solution. 

p. = m./lOOO 
x x 
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5. Compute area density of each element in dry 
mounted standard. 


6. Compute total area density of dry mounted 
standard. 


a 



7. Compute weight faction of each element in dry 
mounted standard. 


w . 

i 


V5X = MIX = a i /a 


8. Compute mass absorption (attenuation) coeffi- 
cients of 1th element for Ka line of k th element 
analyzed. 




c(z.)X 


n (z ) 
k 


(k = 1 , . . . , N) 


c and n are constants obtained from J. W. 
Colby’s set of atomic constants. 


9. Compute mass absorption coefficient of sample 
for each element analyzed. 


x k CT = esc £cr j (p/p i ) k 
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10. Compute absorption factor of standard for 
each element analyzed. 



II. Computation related to the red blood cells 

A. Input data 

-w\ — approximate x^eight fractions of elements in 

dry RBCs (i = 1, . .., L) 

o' — average area density of dry RBC 

N k — accumulated count for k tk element analyzed, 
corrected for background (k = 1, . . . , N) 

B. Computations 

1. Compute approximate average area density of 
each element in dry red blood cells. 

o’. = u.a 

i i 

2. Compute average mass absorption coefficient of 
RBC for each element analyzed. 

X k o = esc ^S a iCv/P i ) k 

3. Compute average absorption factor of dry red 
blood cells for kth element analyzed. 
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4. Compute are^ density o£ k th element 

analyzed using quantitative equation. 



5. Recompute total area density of RBC using new 
area density values of the N analyzed elements. 

N L-N 

° ' XX + ZX ^ 

k=l i = l 

6. Repeat computations (2) through (5) until the 
largest a update is less than some arbitrary 

it 

percent (0.1 percent for example). 

7. Compute the final weight fractions of the 
elements in the RBC. 

a . 

w . - — 

i — 

a 
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3 . USAGE 


3 . 1 Input Description 

3.1.1 Data and Card Specifications 

All integer (I) input is right- justified, alphanumeric (A]) 
input is left- justified, and real number (F) input is 
unrestricted in its specified fields. 


3.1. 1.1 

Card 1 : 

Title and 

Run Information 


FORMAT (A1 

, 14 , 2X , 5A3 

,2X,5A4,2X,6A4) 

Code 


Columns 

Description and Units 

I DATE 

I 

1,10 

Date 

MISS 

I 

11,15 

Mission identification 

MAN 

I 

16,20 

Subject code 

JSET 

I 

21,25 

Sample code 

3.1.1. 2 

Card 2: Program Parameters 

FORMAT (3I5,F5.0) 

Code 

Type 

Columns 

Description and Units 

IPRINT 

I 

1,5 

Namelist; write indicator for 

diagnostic purposes 

(IPRINT 5 0 for "no write") 

IWT 

I 

6,10 

Weight index for subroutine GLHFR2 


preceding 


page BLANK NOT Etu®* 


[0 — weights are input through 
calling arguments as 1 / (sample 
variance). 1 — weights are 
defined as 1 . in the subroutine. 

2 — weights are calculated in the 
subroutine as 1/Y(I)]. 
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3. 1.1. 3 

Card 3: Microprobe 

FORMAT (4F10. 2) 

Parameters 

Code 

Type 

Cc lumns 

Description and Units 

AP 

F 

1,10 

Accelerating potential (Kv) 

BD 

F 

11,20 

Beam diameter CIO " 4 cm) 

CT 

F 

21,30 

Count time (sec) 

DT 

F 

31,40 

Dead time (10 sec) 

3. 1.1. 4 

Card 4: Sample Parameters 

FORMAT (3I5,F5.0) 

Code 

Type 

Columns 

Description and Units 

ITYPE 

I 

1,5 

Sample type (1 — RBC, 2-tissue 
3-gelatin standard) 

NSAMP 

I 

6,10 

Number of samples 

NA 

I 

11,15 

Number elements analyzed 

TBIO 

F 

16,20 

Thickness of samples (10 4 cm) 


J 


3. 1.1. 5 Card 5: Standard Parameters 

FORMAT (3F10. 5,15) 


Code 

Type 

Columns 

Description and Units 

GEL 

F 

1,10 

Gelatin in liter of standard 
solution (grams ' 1 

H20 

F 

11,20 

Fraction of water in stock 
gelatin (0. 5 H20 <_ 1.) 

TSTD 

F 

21,30 

Thickness of standard section 
( 10~ 4 cm) 

NSTD 

I 

31,35 

Number of standard intensity 
counts 
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3.1, 1.6 Card 6: Single Standard Estimated Data 

FORMAT (A2,3X,3F10.0) 


Code 

T/pe 

Columns 

Description and Units 

AE(I) 

A 

1,2 

Chemical symbol for I th analyzed 
element (right justified) 

BG(I) 

F 

6,15 

Statistically estimated back- 
ground for 1th analyzed element 

ESC (I) 

F 

16,25 

Statistically estimated intensity 
count for 1th analyzed element 

SDVESC(I) 

F 

26,35 

Standard deviation of estimated 
intensity count for 1th analyzed 
element . 


Note : Program requires one card 6 for each element analysed. 

These cards are punched out when program runs statis- 
tical analysis of standard data 3 ITYPE = 0 . 


3. 1.1.7 Card 7: Standard Chemical Composition 

FORMAT ( I 1 , 3X , 6 (A2 , IX , F6 . 1 , IX) ) 


Code 

Type 

Columns 

Description and Units 

NAD 

I 

1 

Number elements added to 
standard 

ASE(I) 

A 

5,6 

Chemical symbol of added element 
(right justified) 

AMTASE(I) 

F 

8,13 

Amount added element (meq/&) 

ASE (n) 

A 

lOn -5, 
lOn -4 

nth chemical symbol 
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Code Type Columns Description and Units 

AMTASE(n) F lOn -2, nth amount 

lOn +3 

Note: Program dimension limitation , n _< 6 . 


3.1.1. 8 Card 8: Standard Intensity Counts 

FORMAT ( 1 2 , 2X , 4 (A2 , IX , F6 . 0 , IX , F6 , 0 , IX) ) 


Code 

Type 

Co lumns 

Description and Units 

IPS 

I 

1,2 

Principal standard indicator 
(IPS = 1 designates principal 
standard) 

ADSE(l) 

A 

5,6 

Chemical symbol of first 
analyzed element 

SCT(1,ND) 


8,13 

Intensity count of first 
analyzed element in NDifc/i standard 

SCTSD (1 ,ND) 


15,20 

Standard deviation of estimated 
count 

ADSE(n) 


17n -12, 
17n -11 

For nth element 

SCT (n,ND) 


17n -9, 
17n -8 

For nth element 


Note: Program dimension limitation is n <_ 10 , and card 

space limitation is four elements per card. 
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3. 1.1. 9 

Card 9: Biological 

FORMAT (6 (8X,A2)) 

Elements Analyzed 

Code 

Type 

Columns 

Description and Units 

ABE (1) 

A 

9,10 

Chemical symbol of first 
analyzed element 

ABE (n) 

A 

lOn - 1, 
lOn 

For nth analyzed element 

3.1.1.10 

Card 10: Sample Intensity Counts 

FORMAT (6F10. 0) 

Code 

Type 

Columns 

Description and Units 

CTB [1] 

F 

1,10 

Intensity count of first 
analyzed element 

CTBCn) 

F 

lOn -9, 
lOn 

For nth analyzed element 


Note: Program dimension limitation is n <_ 10 , and card 

space limitation is n <_ 7 for cards 9 and 10. 

3.1.2 Tape Specifications 

Program tape, NASA-MSC Computation and Analysis Division 
No. A06690 has one file containing the program with 
subroutines . 

3 . 2 Program Run Preparation 
3.2.1 Deck Setup 


3-5 


3. 2, 1.1 Punched Program Deck Run 

Running the program from a punched program deck requires the 
following card sequence: 

1. NASA-MSC run card (V RUN) 

2. Message card specifying one tape (VN MSG 1 TAPE) 

3. Microfilm plot control card (V PLT) 

4. Input tape assign card (V ASG A=4483) 

5. BIOMAP control card (V FOR BIOMAP) 

6. BIOMAP program deck 

7. GLHFR2 control card (V FOR GLHFR2) 

8. GLHFR2 subroutine deck 

9. IGIVE control card (V FOR IGIVE) 

10. IGIVE subroutine deck 

11. GROP1 control card (V FOR GR0P1) 

12. GROP1 subroutine deck 

13. PLOT3 control card (V FOR PLOT3) 

14. PLOT3 subroutine deck 

15. HIS1 control card (V FOR HIS1) 

16. HIS1 subroutine deck 

17. PCORRE control card (V FOR PCORRE) 

18. PCORRE subroutine deck 

19. GAMMA control card (V FOR GAMMA) 

20. GAMMA subroutine deck 

21. FACTOR control card (V FOR FACTOR) 

22. FACTOR subroutine deck 
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23. Program execute control card (V XQT BIOMAP) 

24. Data deck 

25. End- of -file control card (V EQF) 

26. Optional post mortem dump control card (VE PMD) 

3. 2. 1.2 Tape Run 

Running the program from a tape requires the following: 

1. NASA-MSC run card (V RUN) 

2. Message card specifying two tapes (VN MS G 2 TAPES) 

3. Microfilm plot control card (V PLT) 

4. Input tape assign card (V ASG E=$STATP) 

5. Program tape assign card (V ASG B= ) 

6. Complex utility program execute control (V XQT CUR) 
and the four following CUR statements 

7 . TRW B 

8. IN B 

9. TRI B 

10. TOC 

11. Control cards for listing or changing the program or 
subroutine statements (V FOR,* NAME, NAME) 

12. Follow with cards described in paragraph 3. 2. 1,1, 
sections (23) , (24) , (25) , and (26) , as for program 
deck run. 
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3.2 r 2 Special Control Cards 

Card 2, the program parameter card, contains (1) IPRINT which 
must be assigned a nonzero integral value to obtain a name- 
list printout for diagnostic purposes, and (2) the weighting 
index IWT to specify the weighting option for subrou- 
tine GLHFR2. 

3.2.3 Special I/O Devices 
None . 

3.2.4 Overlay Structure 
None . 

3 . 3 Output Description 

3.3.1 Printer Output Identification 

A run which analyzes a number of standard intensity counts 
to statistically estimate the intensity count of the primary 
standard will print out a table giving, for each added 
element of each standard, the total milli- equivalent amount 
(residual plus added], the measured intensity count, the 
area density, the weight percent, the absorption factor; 
and, for the primary standard, the estimated intensity count 
with standard deviation, the estimated background count with 
standard deviation, and the estimated intensity count less 
the background. 

For a run analyzing standard samples relative to a primary 
sample a table is printed out giving for each analyzed element 
in each sample the measured intensity count, the absorption 
factor, the absorption correction, the area density, the 


weight percent and the calculated concentration of the 
standard samples. 

For a run analyzing a set of red blood cells relative to a 
primary sample, a table is printed out giving for each 
analyzed element in each red blood cell the measured inten- 
sity count, the area density, and the weight percent; and, 
for the set, the weight percent range. In addition, a 
cumulative plot and a histogram of the weight-percent dis- 
tribution over the set is printed. 

3.3.2 Microfilm Output 

For each analyzed element in a set of standards a microfilm 
plot is made which relates the estimated linear relation 
between the intensity count and the absorption- corrected 
area density for that element. 

3. 4 Execution Characteristics 

3.4.1 Restrictions 

The following analytic restrictions must be observed: 

(1) When analyzing a set of standards to estimate a most 
probable intensity count for the principal standard, 
assign the integer value of one to the program parameter 
IWT if the standard data are single-point microprobe 
intensity counts. If the standard data are average 
intensity counts over samples taken from each standard 
and the sample variances, assign the integer value two 
to IWT. 
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(2) The number of standard data must be greater than two. 

If not, the program will stop in the subroutine GLHFR2 
after writing a diagnostic message noting this error. 

(3) When specific standard or sample data are input 

(card 10) , they are accompanied by the chemical symbol 
of the corresponding elements (card 9) . If these sym- 
bols are not found when checked against the array of 
analyzed element symbols the program stops after writing 
out a diagnostic statement noting this error. 

(4) Principal standard data must be last in the sequence of 
standards . 

3.4.2 Running Time/Lines of Output 

A multiple- execute analysis of sodium and potassium in sets 
of 100, 35, 35, and 35 red blood cells ran 1.5 minutes. A 
combined statistical analysis of five standard averages and 
sodium and potassium analysis of nine standard samples ran 
0.5 minute. Lines of output can be calculated as one line 
per element analyzed in each standard or red blood cell 
analyzed. A run performing all analyses of red cell sets 
outputs cumulative distribution and frequency distribution 
plots over weight percentage, each of which require three 
pages per element analyzed. 

3.4.3 Accuracy /Validity 

Arithmatic operations are single-precision, which is quite 
adequate for the four significant figure iterative criterion 
of the area density computation. This computational preci- 
sion is an order of magnitude better than the precision of 
experimental measurement. 
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4. REFERENCE INFORMATION 


4c 1 Detailed Flow Chart of BIQMAP 

The detailed flow chart of BIOMAP is shown on the following 
pages. 


srecedmg page blank not nnafl 
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4. 2 Symbol Definitions 


Math 

Code 

Type 

Description 

x n 0 

ABC (n) 

F 

Mass absorption coefficient of 
biological sample for nth 
analyzed element. 

V 3) n 

ABCO (n) 

F 

Absorption correction for nth 
analyzed element. 


ABCS (n.) 

F 

Mass absorption coefficient 
of standard for nth analyzed 
element . 


ABE (n) 

I 

Chemical symbol of nth analyzed 
element in biological sample. 

ia 

AC Cn,m) 

F 

Mass absorption coefficient of 
mth element for Ka line n c nth 
analyzed element. 

cr 

ADB 

F 

Area density of biological 
sample . 

a 

n 

ADNB (m) 

F 

Area density of m th element 
in biological sample. 

0 

a 

m 

ADNS Cm) 

F 

Area density of mth element 
in standard. 

a° 

ADS 

F 

Area density of standard. 


ADSE Cn) 

I 

Chemical symbol of nth analyzed 
element in standard. 


AECn) 

I 

Chemical symbol of nth analyzed 
element in single standard. 
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Math 

Code 

Type 

Description 

a 

n 

AFB(n) 

F 

Absorption factor of biological 
sample for nth analyzed element. 

0 

a 

n 

AFS (n) 

F 

Absorption factor of standard 
for nth analyzed element. 


AMTASE (&) 

F 

Amount in millimoles per liter 
of Hth element added to standard 

0 

w 

m 

AMTGEL(m) 

F 

Approximate weight percentage 
of m th element in dry gelatin. 
Literal list in data statement. 

W 

nc 

AMTRBD Cm) 

F 

Approximate weight fraction of 
mth element in normal dry red 
blood cell. Literal list in 
data statement. 

E o 

AP 

F 

Kinetic energy of focused 
electrons given in Kev. Also 
electron gun accelerating 
potential. Not used in calcu- 
lation. 


ARATIO 

F 

Ratio of effective area of 
excitation on the standard to 


that of the biological sample. 
Assigned a value of 1.0 for 
this version of BIOMAP. 
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Math 


2 

ra 


Code 

Type 

Description 

ASE (&) 

I 

Chemical symbol of 1th element 
added to standard. 

BD 

F 

Beam diameter (10 -4 cm) . 

BG(n) 

F 

Background intensity count for 
nth analyzed element in standard. 
Ulote: Concentration intensity 

count for biological samples 
are assumed to be background 
corrected. 

BIOEL(m) 

I 

Chemical symbol of m.th element 
in normal biological cells. 
Literal list in data statement. 

BIOZ(m) 

I 

Atomic number of m th element in 
normal biological cells. Literal 
list in data statement. 

COEF(k) 

F 

k th component of coefficient 


array for given analyzed ele- 
ment returned to BIOMAP from 
GLHFR2. C0EF(1) is the esti- 
mated background intensity count 
(intercept) and C0EF(2) is the 
estimated rate of concentration 
intensity count increase xvith 
absorption corrected area den- 
sity (slope) . 
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Math 


Code 


Type Description 


c 

CON 

F 

Empirical constant C used in 
computing mass absorption 
coefficient (re^ = C> n ) . 
Obtained from atomic constants. 


CONC (n) 

F 

Concentration of nth analyzed 
element in gelatin sample 
(millimoles/liter) . 

Wo 

CR(n) 

F 

Concentration intensity count 
ratio for nth analyzed element 
(sample count/standard count) . 

csctJj 

CSCTH 

F 

Cosecant of 52.5 degrees the 
angle of X-ray emergence. Fixe 
parameter set equal to 1.2605. 


CT 

F 

Counting time for sample and 
standard (seconds) . 

I 

n 

CTB (n) 

F 

Concentration intensity count 
of nth analyzed element in 
biological sample. 


DATE 

A 

Date of computer run request. 


DEL(n) 

F 

Fractional difference between 
area densities of nth analyzed 
element as computed in two 
adjacent iterative cycles. 


PENAD (&) 

F 

Volume density of %th element 
added to standard (grams /cm ) . 


DENS (m) 

F 

Volume density of m th element 
in gelatin (grams/cm ). 
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Math 


Code 


X 




DRYGEL 


DT 

EDGE 


ESC (n) 


EX 


F 


GEL 

HZO 


Description 

Concentration of anhydrous 
gelatin in standard base 
(grams/cm 3 ) . 

Counter dead time (10 seconds) . 
Not used in computation. 

Absorption edge wavelengths X 
used in computation of mass 
absorption coefficients 
(y K = CX n ) . Obtained from 

K 

atomic constants, 

F Statistically estimated concen- 

tration intensity count for nth 
analyzed element of single 
standard. 

F Exponent n of absorption wave- 
length X used in computation 
of mass absorption coefficients 
(y :C = CX 1 " 1 ) . Obtained from 

fC 

atomic constants. 

F An array of weight percentage 

frequencies falling in each 
interval or group. Returned 
from statistical subroutines 
I GIVE and HIS1. 

F Concentration of suock gelatin 

in gelatin standard (grams/liter). 

F Fraction of water in stock 

gelatin. 


Iffi! 

F 

F 

F 
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Math 

Code 

IXE£ 

Description 


I COUNT 

I 

Iterative cycle counter for 
sample element area density 
computation. 


ICS 

I 

Array of chemical symbols for 
the chemical elements. 


INDX 

I 

Holder for sample number of 
principal standard. Always 
last standard in sequence. 


IPS 

I 

Flag identifying principal 
standard. Always last standard 
in sequence. 


ITITX 

I and 

Data statement arrays containing 


ITITY 

I 

abscissa and ordinate titles for 
microfilm plot of concentration 
intensity count vs absorption 
corrected area density. 


ITIT1 

I and 

Data statements containing 


ITIT2 

I 

histogram of print plot headings. 


I TYPE 

I 

Type designator for sample to be 
analysed (RBC-1, TISSUE- 2, 
GELATIN- 3) . 


MDIM 

I 

Dimension of concentration inten- 
sity count and absorption cor- 
rected area density array input 
to subroutine GLHFR2, 


NA 

I 

Number elements analyzed. 


NAD 

I 

Number elements added to standard. 
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Math 


Code 


Type Description 


NAME 

A 

Name of run requester. 

NBE 

I 

Number of biological elements. 

NG 

I 

Number of weight percentage 
groups. Fixed in program. 

NPROB 

I 

Problem number in title and 
run information (Card 1) . 

NSAMP 

I 

Number of biological samples 
to be analyzed. 

NSE 

I 

Number of elements in standard. 

NSTD 

I 

Number sets of standard concen- 
tration intensity count data 
used in statistical estimates. 

NZ 

I 

Atomic number symbol used in 
computation of mass absorption 
coefficient. 

NZA(&) 

I 

Atomic number of J ith chemical 
element added to standard. 

NZAE(n) 

I 

Atomic number of nth analyzed 
element in standard. 

NZSE (m) 

I 

Atomic number of m th chemical 
element in standard. 

P(i.j) 

F 

j th atomic constant of i th 
atomic number chemical element 
(i = 1,100 and j = 1,36). 

PABC 

F 

Partial absorption coefficient 
of standard or sample. 
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Math 


Code 


POE 


Q 

R 


RHOB2(n) 


SCT (n) 


SDEV 

SDVERR(n) 


SDVESC(n) 


SE 

SIGBl(n) 
SIGB2 (n) 


Type Description 

F Power of e (base of natural 

log) ; used in computing absorp- 
tion factors . 

F and Arrays containing coordinates 
F of endpoints of microfilm 

plotted line fitted to concen- 
tration intensity count vs 
absorption corrected area 
density data. 

F Computed volume density of nth 
analyzed element in gelatin 
sample . 

F Concentration intensity count 

of nth analyzed element in 
standard . 

F Standard deviation of RBC popu- 

lation distribution. 

F Estimated standard deviation of 

error in intensity count measure- 
ments for nth analyzed element. 

F Standard deviation of error in 

estimated concentration inten- 
sity count for nth analyzed 
element . 

A Chemical symbol array for ele- 

ments in gelatin. 

F and Area density of nth analyzed 
F element in sample as computed 
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Math 


Code 


Type 


o 

n 


X 


K 


SIGS(n) 

STEP 


SUM1 

SUM2 

sx(ij) 


TAG 

TBIO 

TITLE 

TSTD 

VAR 


VARCOV(i,j) 

VARERR 

WAVE 


F and 
F 


F 


A 

F 

A 

F 

F 


F 

F 

F 


Description 

in current and previous iterative 
cycles, respectively. 

Area density of nth analyzed 
element in standard. 

Width of weight percentage 
groups in cumulative and fre- 
quency distribution plots. 

Cumulative sums used in comput- 
ing RBC distribution statistics. 

Absorption corrected area den- 
sity of j th element in ith 
standard. 

Problem type identification. 

Thickness of biological sample 
CIO -4 cm) . 

Descriptive title of run. 

Thickness of gelatin standard 
(10 -4 cm) . 

Variance of weight percentage 
distribution of analyzed ele- 
ment in RBC sample. 

ij element of covariance matrix 
returned from GLHFR2. 

Estimated measurement error 
variance returned from GLHFR2. 

Wavelength of Ka X-ray photon 
emitted by analyzed element. 
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Math 

Code 

Type 

Description 

w 

m 

WFBE (mj 

F 

Weight fraction of m th element 
in biological sample. 


WM 

F 

Array of statistical weights 
used in GLHFR2. 


WTPMX(i,s) 

F 

Weight percentage of j th ana- 
lyzed element in ith RBC . 


XCi) 

F 

Weight percentage of specified 
element in i th RBC. 


XBAR 

F 

Mean weight percentage of speci- 
fied element in RBC set. 


XL ,XR 

F 

Left and right marginal values 
of absorption corrected area 
density in microfilm plot. 


XM 

F 

Array of absorption corrected 
area densities. 


XMAX 

XMIN 

F 

Maximum and minimum weight per- 
centage of specified analyzed 
element in RBC sample. 


XMX 

F 

Maximum absorption corrected 
area density of specified ana- 
lyzed element in set of standards 


YB ,YT 

F 

Top and bottom marginal values 
of intensity count for micro- 
film plot. 


YM 

F 

Array of intensity counts for 
specified analyzed element 
over set of standards. 


YMX 

F 

Maximum intensity count of speci- 
fied analyzed element in set of 
standards . 
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4 . 3 Subprogram Documentation 


4.3.1 STAT-CAT Statistical Subroutine 

Documentation of the following subroutines can be found in 
the NASA-MSC Computation and Analysis Division Statistical 
Catalogue: XGIVE, GR0T1, TL0T3, HIS1. 

4.3.2 Subroutine GLHFR2 
Documentation follows. 
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SUBROUTINE GLHFR2 


IDENTIFICATION 

Name/Title 

Author/Date 
Documented by/Date 
Organization/ Ins tall at ion 
Machine Identification 
Source Language 


— GLHFR2 (General Linear 
Hypothesis Statistical Model 
of Full Rank Two) 

— B. F. Edwards, March 1971 

— B. F. Edwards, March 1971 

- LEC/NASA-MSC 

- UNIVAC 1108 

- FORTRAN V 


PURPOSE 

The purpose of this subroutine is to compute best linear 
unbiased estimates (BLUE) of parameters associated with a 
linear statistical model of full rank two. This model has 
the form y = Xg + e where Y is a random n-vector of 
weighted observations (measurements) , X is an nx2-matrix 
of weighted known fixed independent quantities (experimental 
parameters) , fj is a 2-vector of unknown coefficients 
(intercept and slope) , and e is an unknown random n-vector 
of measurement errors which satisfies either Case (1) : 
e is distributed N(Q,o 2 I), with o 2 unknown, or Case (2): 
e is a random vector such that E(e) = 0 and 
cov(e) = E(ee') = a 2 I , with a 2 unknown. 


3/71 


GLHFR2-1 
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USAGE 


• Calling Sequence 

CALL GLHFR2 (N ,X,Y,W,M.EB , EVAR, VAREB) 


Math 

Code 

Type 

Description 

n 

N 

I 

Number of observations. 

tx.} 

X 

F 

n- array of known experimental 
constants . 

{yp 

Y 

F 

n* array of observations. 

{Wi } 

W 

F 

n-array of weighting coefficients. 


M 

I 

Calling program dimension of the 
X, Y, and W arrays. 

(3 

EB 

F 

2-array of linear coefficient estimate. 

A 2 
0 

EVAR 

F 

Estimate of equivalent homegeneous 
error variance. 

/s 

covCB) 

VAREB 

F 

4-array of variance- covariance matrix 

A 

elements for cov 8 in the order: 11, 


12 , 21 , 22 . 

« Data In/Out 

Data is input or output through argument lists and the 
printed namelist HATS which contains N and EVAR and 
the elements of EB and VAREB. 

• Error Messages 

An error message is printed out and the program stops 
if N is not greater than 2. 

3/71 GLHFR2- 2 
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• Storage 

Code: 252 octal 

Data: 143 octal CM = 25) 


METHOD 

I£ a quantity to be measured (y) is known to depend 
linearly on a quantity (x) over a given range (x f , x") , 
the functional relationship between y and x is given 
by 

y = a + bx ; (x 1 < x < x") (1) 


where a and b are constant parameters. If a pair of 
points (x , y ) and C x 2 ’ l^) are known which satisfy the 
relation (1) , then a and b can be computed. 


If, however, for fixed x , say x. , one successively 
attempts to measure y(x i ) but obtains fluctuating values 
y where the subscript j denotes the j th measurement 
at x i , then y can be equated to the true value 
y(x^) plus an error , 



pCXf) 


+ 



or 


y 


ij 


a + bx. + 

X 



( 2 ) 
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Solving Equation (2) for the error term yields 


e., = y.. -a-bx. 

13 7 in 1 


C3) 


I£ (A), the e.. are normally distributed with mean 
1 3 o 2 

0 and variance cn , that is e_ — 14(0,^}, or (B) , 

the expectation of is 0 and the expectation of 

(e . . ) 2 is a 2 ; and, furthermore, for either (A) or (B) 

the erf are equal for every i , then point and interval 

estimates of the parameters a and b and y (x) can be 

made provided that the total number of measurement is 

greater than two and i > 1 . 


If the o 2 are not equal, e„ in Equation (3) may be 
transformed to e* . ~ N(0,1) so that all a 2 are equal 
to unity. This transformation consists of multiplying 
(3) by 1/a to obtain 



y . ./a . - (a + bx . )/a . 
7 x;j x v 1 x 


In general the ck are not known and estimates of o ± , 
s. must be obtained from independent measurements or as 
functions of y. 

e* . - y. ,/s. - (a + bx.)/s. (4) 

X3 / 2.3 1. x x 

with e* . distributed with variance a* 2 only approxi- 
13 

mately equal to one. 
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The n measurements (4) can be combined into the matrix 
equation 


e = Y - X3 (5) 

where e and Y are n-vectors; 

X is a nx2-matrix with row elements f/w7 , x . /wT) where 
w i , the weighting factor, is 1/s^; and 3 is the 2-vector 
of parameters a and b . Estimates of a and b can be 
obtained from estimators derived, in case (A) by the maximum- 
likelyhood method or, in case (B) , by the least- squares 
method. In the least-squares method 

e* 2 = e'e = (Y - XB) ' CY - XB) (6) 

Differentiating with respect to & and setting to zero 
yields 

B(e , e)/33 = 2X’Y - 2X'XB = 0 

A 

and the unbiased estimate of 3 , 3 is the matrix solution 

3 - (X , X)“ 1 X , Y (7) 

An estimate of a , the expectation of (e^) following 
from (6), is a* 2 ~ [Y - X^} ' (Y - X3)/n . Slight modifica- 
tion produces an unbiased estimate: 

a* 2 = (Y - Xer CY - X3)/(n - 2) 

or 

(Y T Y - 3 ' X ' Y) / (n - 2) (8) 

3/71 GLHFR2- 5 
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Expansion of (7) yields: 


with 


a = 


b = 


D = 


s i = (EVi Ev 2 ‘ E'Vi I>i x i y i)/ D 
e 2 = (2>i Ev^i - E w i x i E'^yJ/o 

E«i E w i x i - (E'Vi ) 2 


Expansion of (8) yields: 
"*2 


■ (Ev^ ‘ a E» 1 y i - bE w i x t y< ) 

A 

The variance-covariance matrix of 3 ; s given by 


and yields 


(9) 


(10) 


A 

cov(3) 

“ 

(X'X)" 1 ^* 2 = (X’X) 


"2 

O' 

a 

A * 2 
a 

2>i x i /D 


2 

ab 

"2 

C ba 

= -a* 2 2 w i x i/D 


11 

w a 

< to 

;* 2 

E«i/D 

(11) 


In this subroutine the n-arrays of measurements Y , 
known fixed constants X , and precalculated weighting 
coefficients W are input through the calling arguments 
and the estimates in equations (9) , (10) , and (11) are 
calculated and output, also through the calling arguments. 

3/71 GLHFR2- 6 
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• Symbol Definition 


Math 

Code 


Definition 

D" 1 

DETINV 

F 

Reciprocal of the determinant 
D = | X T X | = ] S | . 

(X’X)” 1 

SINV11, 

SINV12, 

SINV22 

F 

The 11, 12, and 22 elements of 
the (X'X) -1 or S -1 matrix. 

w . 

i 

ww 

F 

Working symbol for the weighting 
coefficient w. . 

X 

X . 
X 

XX 

F 

Working symbol for the fixed 
parameter x . 


YY 

F 

Working symbol for the measure- 
ment y. . 

w.x, 

X X 

wx 

F 

Product symbol. 

2 

W.X . 

X X 

wxx 

F 

Product symbol. 

w . x . y . 

X x' X 

WXY 

F 

Product symbol. 

w.y , 

x' X 

wy 

F 

Product symbol. 

2 

i^r.y . 

X X 

WYY 

F 

Product symbol. 

Ew. 

SUM 

F 

Sum of the w. . 

X 

2 w i x i 

SUMWX 

F 

Sum of the products w.x., 

XX 

\ . 2 
Vw . x . 

x - r X X 

SUMWXX 

F 

2 

Sum of the products 

T) w . x . y . 

1 f X x'x 

SUMWX Y 

F 

Sum of the products of w 
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Math 

Code 

Type 

Definition 

Dvi 

SUMWYY 

F 

2 

Sum of the products w i y i « 

5> 2 i 

SUMWW 

F 

Sut 1 of the producti w? . 

X'Y 

XTY1 , 
XTY2 

F 

Components of the vector (X'Y). 


• References 

F. A. Graybill, An Introduction to Linear Statistical 
Models , Vol. 1, McGraw-Hill Boo? mpany, Im.. , 

New York (1961) , Chapter 6. 

C. A. Bennett and N. L. Franklin, Statistical Analyses 
vn Chemistry and the Chemical Industry 3 John Wiley and 
Sons, New York (1954), Section 6.27. 

RESTRICTIONS 

• Analytic 

The number of observations (measurements) must be 
greater than two. 

• Hardware 
None. 

• Operational 
None. 
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ACCURACY 


The subprogram uses single-precision accuracy. There is 
no loss of significance in the output. 

VALIDATION PLAN 

Mathematical equations were checked out by showing, when 
the weighting coefficients were equal to unity, that they 
reduced to the least-squares line-fitting formula. Opera- 
tional checks were made by showing that the straight line 
estimates correctly fit the data sets. 

CODING INFORMATION 

• Special Program Constants 
None. 

• Timing 

The subroutine has not been timed. 

DETAILED FLOW CHART 

• See the following page. 
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MSC Form 887A (Mar 65) 










SOURCE LISTING 

See the following pages. 


3/71 GLHFR2-11 


4-27 
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■ lo.’sinu is? 


»■ FUR,* — - ■ 

UN 1 V AC HUB FORTRAN V LEVEL 2206 QU2Ha - 1 EXEC a LEVEL E12O10C1O9A) 

THIS CuHPl L ATlDN AtS DONE ON |V MOV 7l AT lOi^lill 


subroutine (iLhfr2 entry point 000291 


STORAGE U5Ep! eoOEU-l 00030&I DAT A ( U I 

'SS,: . 

EXTERNAL REFEg^CES (BLOCK, NAME) 


0Uu3 NflOuS 
OUQft N1Q2& - 

00 Ufa NftNL* 
OUB* hSTOP* 
UU07 N E R K 3 S 


000 l 5Qi 


blank COMMON ( 2 ) 000000 


STORAGE ASSIGNMENT {BLOCS, TYPE, RELATIVE LO 


UU0U 
□ Uti 1 
OUti-U- 

- Lt— 

000065 

00021ft 

rruni’i! 1 . 

|E 

13L 

. ry r T 1 , [ . 

OOUl 
OQU 1 

Uruin 


OOOUSO 
0001)45 
.. n n n ■ 1 i "t 

10L 
lftHG 
HATS 
5 1 N V22 

UU0U 

rv 

R 

U UtltT » / 
00002 | 

yt 1 i fiT 
SlNV 12 

OOUO 

R 

uCfTJU J J 

0D0022 

Q0L|0 

K 

□0000b 

SUMTiXY 

UOUO 

R 

OUQU03 

SUMWT 

uuuu 

H 

OUOu 1 ft 

nxx 

OOUO 

R 

Q00U15 

HXt 

Ub(jU 

ft 

00002* 

X T Y 2 

uquu 

R 

OOOU07 

XX 


CATION, NAME) 

□001 000035 1176 0001 000205 ) 2u 

0000 OOOIOS 2F 0000 O0007ft JF 

0000 I 00000a— 1 OOQO oo&l-z-l— fNJPS- 

0000 R OOoOOl SUHW ODQO R 00n002 5UMWX 

onoo R 000006 SUM'RTy 0000 R 00001 1 fjW 

onoo R 000013 AT 0000 H DOOoU Y(YY 

0000 R nOOOlD TT 


0001 0000*15 1276 

ODD! 0000*12 8l 
0000 K 000020 SInV-1 I - - 
onoo H QO 0 OOH 5UHWXX 
0000 R OOOO 1 2 WX 
0000 « OOon23 XTTl 


00 IU1 

1 * 



SUBROUTINE GLHVR2 (N,X ,Y,«V,M,lft,EBiEVAR , VAREB 1 

OOIUS 

2 * 

C 

« 

ft * number uf observations 

00 1 u 1 

3* 

c 

ft 

a = array of known independent variable values 

00101 

ft* 

c- 

■ 

t- * — array- — of observed dependent variable- VA-toes 

0 0)01 

b« 

c 

ft 

A x ARRAY of weight values 

ODtUl 

t> * 

c 

ft 

H X DIMENSION OF INPUT ARRAYS1 K,Y,*. 

UU1U1 

7 • 

c 

9 

EB * ESTIMATES OF LINEAR COEFFICIENTS 

QO 1 0 1 

B* 

c 

ft 

EVAR * ESTIMATE OF ERROR VARIANCE 

00101 

9 * 

c 

ft 

VAREb * elements of variance HATRIX for estimated LIN. COEFFICIENTS 

uotoi 

ItF* 

-c 

— - 



ootui 

u* 

c 


Irt =* 0 AEJomTS are INPUT THROUGH CALLING ARGUMENTS U/5AMPLE VAR) 

uuiul 

u* 

c 


IA ■ l WEIGHTS ARE DEFINED AS 1. IN GLMFR2 

OOIUI 

13* 

L 


Hr * 2 WEIGHTS ARf CALCULATED IN GLHFR2 AS 1,/YIJ) (IFPOISSOm VAR) 

00IU1 

H* 

c 



0U103 

IS* 



DIMENSION X (Ml , Y(nl ,JY(H) ,E8(2) , ViREB (2 ,2) 

OU | OH 

16* 



F&T/f+ATS/W.EtH 1 ) i£B*2) , E V A R ,VAREBt ! , 1 ) T V AHE 0l H 24-rV ARE9 1 2 , 1 1 T 

OUlUft 

l 7 • 



l V ft REB ( 2 i 2 1 

OU 105 

lb* 



tiRlTEi6»l * 

00107 

19* 



l FORMAT! IhU,///, |0 X,*SUBROUTINE GlHF R 2 ' , / / / 1 

UO t 1U 

2U* 



IF tN.LE.2l GO TO | j 

001 12 

Zl« 



IFlU.EO.l) GO TO S 

001 1H 

t’f 



1 F T Ff'it, r , 1 1 GO To 10 

00114 

23* 



On 7 i = i , n 

0 012 1 




th 1 1 = 1 ,/rtn 

001 22 

2 S* 



/ CunT ) DUE 




6Z-t> , 

ST-ZHfjmD ' . T i/c 


00+21 
0012b 
OP 1 26 
UPU1 
OP l 22 

□ 0 131 

-00+35-— 

00136 

□ 0137 

001 <10 
OP 111 
00112 
OOH-2 

□ 0113 
00116 
00117 

00150 

□ 0151 
-UQ152 - 

00153 

00151 

00155 

00156 

00157 
00160- - 
00161 
00162 

UO 163 
00161 
00 161 
00166 

00 167 

00170 
U0171 

00172 

001 7J 
00173- 

00173 

00171 
00175 
0U176 
00177 
U02U0 
002U } 
002U2 
0P2U2 
00203 
00206 
1)020 7 

002 I 0 
00212 
00213 
00211 
002 1 5 
1)0 2 1 t 
DU22U 
00221 


26* 

27* 

2a* 

29* 

30* 

31* 

32* 

33* 

31* 

3b* 

36* 

37* 


50— T-0— J-U 

a CONTINUE 
DO 9 1 - l i N 
6(1) * 1 . 

9 CONTINUE 
lO CONTINUE 


-J P B 5 - 


39* 
ID* 
11 * 
12 * 
13* 
-H*— 
15* 
16* 
17* 
Id* 
19* 
50*- 
51* 
52* 
53* 
51* 
55* 
Lrirl 
57 * 
Sti* 
59* 
60* 

6 1 » 
-62* - 
63* 
61* 
65 * 
66 * 
67* 
66 * 

6 V* 
70* 
71* 
72* 
73* 
71* 
75 * 
76* 
77* 
7d* 
79* 
dtt* 

til* 

62* 


Sun«x 

sti^n- 

sxJhftxx * 
— SUttftXY * 
SUMftYY * 

DO 1 1 l * 
XX * Kill 
YY * Y t 1 1 
ft* ■ All) 

WX ** 6ft • 


ar~ ti i ~ 

* o* 

* 0. 
» 0 , 
* 0 . 
« 0, 


1 |H 


-ft-f- 


-ftW— *• 


XX 
Y-Y— 

ft Xx a ftx • XX 
TiXY * ftX • YY 
6 Y Y * ft Y • tY 
SOhft * 5u»ft + 

5 U n ft X ■ SOHftX 


ftW 

* WA 
-ft Y 

5UHWXX + ftXX 
SUHftXY + ft X Y 
SUMWTY + ft Y Y 


SUHftXX =• 

SOHrtXY " 

5 U ft ft Y Y * 
ll CONTINUE 

> form s-invekse matrix is 


— DE T I N V" 
5 I N V I 1 
5 l M V I 2 
S1NV22 


* X*XI A NO X 1 Y VECTOR 


-+T/+5UHft*5uMWXX - SUttftX.gUMftX) 
SUHftXX • DET1NV 
“SUKftX ■ DETIMV 
SUtlft * DET1NV 


XTY 1 -m SuMnY 
X T Y 2 a SUNnXY 

VECTOR, ER-ROR- VJHH AHCE* -AND V ftR-|-*ne£-“t u V-AFH AF*CE— Mft-TfH-X 

for beta vector 

EttH) » 5 I NV 1 1 • X T T | ♦ S|NY12*XTY2 

50(2) » SlNvl2*XTYl * S I NV 22 • XT Y 2 

EVftR * (SOtrnYY - (EBU)*XTYI + ER ( 2 ) *XT Y 2 )) /FLO AT « N - 2) 

VAREB t 1 i 1 I * F V Ar • S INV I 1 


V*frE*|- H - r2 1 * -fa VAR * g 1 NV 1 2 ■ - 

VAREB (2*1 1 * V A REO ( I ,2) 

V ARtB ( 2 j 2 ) * EVAR * S1NV22 
, PRINT INPUT AND OUTPUT DAT* 
ftHiTp(6,>|AT5) 

GO To 13 

12 eOttTENOE — — 

«*RlTE(6»3» 

3 F*)H 1 *T( Ihl ,/// I lux , *ERROR - N MUST DE GREATER THAN 2 ») 
STOP 

t 3 CONTINUE 
AK 1 TE < 6 , 2 f 

2 FOrtHft TI two | / / » 1 Ox I • RF TURN T 0 CALLING ROOT t NE *,///// ) 
RETURN 

end 



3/71 ! , j I GLHfRZ-14 


NO b I A fiHOS I4CS 


^ auHF«2 

SYMBOLIC 

if NOV 71 

10:20133 

0 

UIM12HS 

14 

82 

(DFLFTEDl 

*”■ OLHFK2 COOt 

relocatable 

If NOV 71 

10:20:33 

1 

01&HH2M 

1 8 

l 

IoFLETED) 





0 

0 S 3 N f t 

If 

2fi 




4.4 Program Listings 


The program listings are shown on the following pages. 
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Z£-> fr 


_ J_9. „Nfc¥— 7 l- . 


-■ — ^ QHaP i B 1 O H AP — — " ....^ - - — * 

UN t V AC 1100 F OH TRAN V LEVEL Z2U6 0021A - [EXtCB CEVEl E 1 2fl 1 OOQVA ) 

This CuHP t (. A T I ON-iVAS DONE -ON- 1? NOV 7[ AT 


laiini' 


HAIn PROGRAM 

STORAGE USED! CODE <11 0027 | U ! DAT A ( U ) 015011J BLANK COHMOHtZ) 1)00000 

common blocks: 

U0U3 P Q07DZr, 


tXTE.HNAL REFERENCES (^LOCK, NAME) 
DOOl GLMFR2 

OUOS OUIKMP 

■ 00&6— i itt I'KMt 

OU(J 7 MiS| 

00{0 TSlVE 
00 J l CORAN 
00 J 2 PCORRe 

0U13 URDUS 

*■ — 00-{N frKH-S — * ' — 

0015 N E Q2i 
001*. NrtDUS 
00 J 7 EXP 

0020 NVrrtLS 

0021 NEXP65 

— ■ ■ 0022 StiRT— ‘ 

0U2J NriUCj 
0U21 NSTOPl 


STUKAGe ASSIGNMENT (QlOCK, T TPE t RELATIVE LCCATIONj NAMEl 


00 0 u 

0H145 

IU1F 

0000 

011157 

1 02F 

oncu 

001522 

10271 

□ ODO 

011101 

l 03F 

□ 000 

01 il 76 

1 OIF 

OOtJl 

0glb3? 

| OR 1 G 

0001 

□0 1 63a 

1075& 

0001 

000527 

1 1L 

0001 

001636 

1103G 

□ □Ol 

001171 

1 1 IL 

UUyl 

□01*57 

11170 

0001 

00QS66 

1 1 2L 

000 1 

001666 

1 I23G 

000 1 

001031 

1 I 3L 

ODDI 

001702 

I 1 3*G 

UUol 

001 71 1 

1 use 

0001 

001 726 

1 I5 ? G 

000 1 

001573 

1 I6L 

000 I 

□01737 

1 165g 

0001 

00 1 762 

1177g 

UUQ 1 

U0051H 

UL 

0001 

001 765 

1 ZO^G 

onoi 

002003 

1ZI3G 

000 t 

□02031 

1 2 25g 

oooi 

002O37 

1230g 

0001 - 

- -01120-72 ~r g M6tr 

ooar 

T^^roS 

125^0 

nnot 

0021 1 1 

1 26l<»- 

— 000 1 

0W2t1u 

j?75fi - 

— nnnt- 

ooo^ 63 

13L 

auui 

002152 

1307G 

OOUi 

0\ >70 

I3 20& 

0001 

002176 

13265 

0001 

002250 

1 356g 

0001 

002266 

1 37 1 Q 

oo 0 i 

000550 

lH. 

0001 

00. .'24 

|1'5G 

onoi 

OODOlO 

I12G 

000 1 

002355 

1 130G 

nooi 

002377 

)1l5g 

UUul 

002107 

tlb&G 

0001 

002 \ '7 

1 « 7 0G 

nooi 

002517 

1 & 1 1 e 

OQQl 

002550 

|617g 

0001 

002615 

1 53Sq 

OUgl 

000U26 

1 5HG 

0 DO 1 

002633 

1 51 1 G 

anoo 

nllQS 2 

2F 

aooo 

0110&5 

202F 

ODOO 

0 1 1016 

203F 

OOgO 

OHOiU 

2Q1F 

□ DUG 

0 1 10^6 

205F 

nooo 

0(1103 

206F 

0000 

011032 

20 7f 

oooo 

0tia20 

20BF 

OOuO 

OtN|j2!r-2HlF — 

— — oooo - 

OtlZH 

2 HP 

nnoo 

Ol 123 ! 

2F2F 

0000 

0 11211 

2j3f 

oooo 

□11313 

2|5F 

OOgO 

01133g 

2 1 6F 

uouo 

0 1 HI 1 7 

221F 

0000 

011127 

222 f 

OOoO 

nuns 

223F 

aonn 

0 1 151 7 

221F 

OUgO 

Q l 466 t 

225F 

OOUO 

0 1 467 3 

226F 

0000 

011726 

227f 

0001 

000660 

2 3U 

oooo 

Q1 1531 

23 1 F 

OUgO 

014606 

2 J2F 

oooo 

0 1 H 6 1 6 

233F 

nnoo 

011113 

21 If 

0001 

000651 

2Sl 

onni 

000)62 

25 36 

cool 

000177 

2*5G 

OOUO 

01r377 

3F 

nnoo 

011075 

305F 

oool 

□00215 

31 1L 

oool 

□00306 

312L 

O0 0 l 

000166 

32l 

0001 

D00170 

321L 

nnot 

000315 

3Z2G 

oral 

000361 

3?7& 

0001 

oonina 

31QP 

UOgl 

0003 Ts 

3Sl 

OeO 1 

□00122 

3535 

0001 

000111 

365S 

0001 

000151 

3750 

0001 

002662 

3Sl 

UUul 

U l) 2 6 6 1 

3VL 

OOUO 

U I H u 3 6 

IF 

onoi 

000501 

115G 

0001 

0005G1 

1200 

0001 

0DP762 

11L 

UuUl 

000535 

142& 

0001 

ooqsid 

1155 

onoi 

001022 

1&U 

0001 

000777 

17L 

0001 

00n573 

1706 

00 g fc 

001003 

iBL 

tlDOl 

□ Oiuto 

1?L 

noon 

011062 

5F 

□ 00 l 

0 U 1 0 1 3 

&«L 

□noi 

00n607 

5016 


-33 


UOUL 


- 600 + 01 - 

6-Hfi — 

OUU 1 

— 

-GOB* 3 2 

51+5 

0001 


000*72 

*3 66 

-OO0F 


000763-51*6 

0001 


000711 

5556 

UDy 1 


000733 

562 G 

Oou 1 


000713 

566G 

nnoi 


001120 

S7L 

OOUO 


011115 

SRF 

OOOO 


Dll01‘ 

6F 

ouoi 


001*62 

b 1 L 

0001 


001017 

6265 

000 I 


001052 

6 3 26 

0001 


001720 

61L 

0001 


0010711 

612G 

UtlUl 


001 7 2 1 

b5L 

000 1 


001106 

6S1© 

0001 


001111 

6&1G 

0001 


001130 

66 7G 

0001 


0 0 1 7 G 6 

*7l 

Ouuo 


011030 

7F 

0001 


DO 1 7 97 

70L 

onol 


001 ISO 

7Q3g 

0001 


0U12H 

72 l G 

0001 


001235 

721G 

aoui 


OOZUs 

73L 

0001 


0013Q0 

7 37G 

0001 


001361 

767g 

0001 


0U2D16 

77L 

oooo 


01 ml ‘3 

bf 

■ 6ut) + 


■floznH 

*B'lX 

otio r 


"013 21" 6 1 

8*t,~ 

nno t 


0023 13- 

“9 BE 

— O&01- 

/ Y7L 

- 0001 


00270+ 

•799L" 

OUoO 

ft 

01 1363 

ABC' 

OoOQ 

R 

U 1 1 107 

ABCO 

nooo 

R 

01125s 

ABC 5 

0000 

i 

006061 

ARE 

O0OQ 

ft 

010663 

A ' 

0U()0 

ft 

013627 

ado 

OoOO 

R 

Ol 1 203 

admb 

onoo 

R 

011551 

AONS 

oooo 

R 

013531 

ADS 

onoo 

I 

000032 

At 

uuuo 

ft 

01 137b 

afb 

Oquo 

R 

011301 

AFS 

□ooo 

ft 

0 t 1 5D3 

AMTaSE 

oooo 

R 

011157 

amtgel 

OOOO 

K 

01 1171 

AHYRbC 

OOuU 

ft 

000110 

AM I 0 

ttouo 

I 

0QDU76 

AWSS 

nnoo 

R 

013517 

AP 

oooo 

R 

013530 

arat 1 0 

noon 

I 

00Dn21 

A S'_ 

OUUU 

. . III 1. JL 

ft 

013S&0 

BO 

uouu 


013/31 

BOAT A 

onoo 

R 

011313 

■ i 1 r._ . 

BG 

oooo 

ft fi -. n 

1 

OUODOO 

BIOEL 

naoo 

1 

D 

oooo i z 

BIOZ 

UUUv~ 

-ft~ 

0 1 Jl* J~ 

ce 

thjt)s- 

ir 

-©+0*2+ 

cocf ■ 

onoo - 

ft 

ijr«37o cro rr — 

truw U” 

K~ 

-fl 1+53-?— 

-CON* - 

oooo 

M 

63 3*0 1 

CP 

OOuD 

fi 

0 1 L35l 

CR 

oqou 

R 

01351 1 

CSCTH 

onoo 

R 

D13S5I 

CT 

onoo 

R 

01 1325 

CTB 

oooo 

R 

013+11 

cv 

oooo 

ft 

01 JnZl 

C 1 

Oouo 

R 

0 l 362 & 

C12 

□ riOO 

R 

013526 

C2 

oooo 

ft 

01 1133 

DEL 

onon 

R 

Dl 1517 

DENAp 

ao u o 

ft 

011113 

DENS 

ouua 


013622 

OIST 

onoo 

R 

0 I 3566 

DRY6EL 

□ ooo 

ft 

0 1 3552 

OT 

oooo 

R 

013575 

EDGE 

UOuU 

ft 

011337 

ESC 

OQOO 

R 

007130 

ESCMA 

onoo 

R 

01357/ 

E a 

O0DO 

R 

012521 

F 

onoo 

R 

013557 

GEL 

OUuU 


0 1 3 Z02 

GUM I 0 

Onoo 

R 

013660 

H20 

□ ODD 

1 

013523 

1 

OOCD 

1 

013533 

ICOOnT 

0003 

I 

OOOlU 

ICS 

0ut)0— 1 — 0 1 36 Zt — rO-A T£— 

01300-1—0+3*02 

+ N0* 

onoo- i- 

013&15 

J PH- PUp- 

0006- I 

6-13561- 

JP5 

— — OOQ0 

1 

010*33- 

17}+* 

(JUDO 

1 

0IU64 1 

I T 1TY 

UUUO 

1 

013563 

ITYPE 

QQOO 

I 

0 135lb 

llT 

oooo 

1 

□13565 

J 

OODO 

I 

013621 

JA 

UUuO 

1 

013611 

JSET 

unou 

I 

013570 

K 

noao 

I 

0 1 357 1 

K 

ooou 

1 

013513 

HAN 

OOOO 

1 

013603 

hdjk 

ouou 

I 

013512 

HISS 

OOuu 

1 

013667 

ft 

oooo 

1 

013555 

N A 

oooo 

I 

013563 

NAD 

OOOO 

1 

013611 

NDE 

dutfO 

1 

013662 

NO 

OflUD 

1 

0 1 3 6 i 7 

NG 

0000 

1 

013613 

NS 

oooo 

1 

013551 

NSAHP 

OOOO 

1 

0)3571 

NSE 

00y0 

1 

013&32 

ftSTO 

Oouo 

I 

013573 

NZ 

DnOO 

1 

011511 

NZ A 

onoo 

1 

01 1213 

NZ At 

oooo 

1 

01| 123 

NZSE 

OUtji 

ft-oeottotr 

-P 

&W9-H-013WC0 

pab-c 

- npoo 

R 

013*20 

-PCOR 

66+2- 

-R- 

-ODOOOO-Pe-OftRE 

- 6000 

K 

013*01 

POE 

ouou 

ft 

0 10b 3 1 

ft 

UQUO 

R 

U 10627 

r 

onoo 

R 

01 1 12 I 

hhobz 

0000 

R 

010111 

5CT 

oooo 

R 

011615 

SCTSO 

OOuU 


013616 

soaTa 

oouo 

R 

01 35i0 

soev 

onoo 

R 

Dl 1603 

sove r r 

□ ooo 

ft 

0HS71 

snvEst 

oooo 

1 

0OGQ1+ 

5E 

auou 

ft 

Oil lib 

S16B1 

uoou 

R 

01 1625 

SIGB2 

oooo 

R 

011267 

S IGS 

D0aa 


013771 

5TAT 

oooo 

R 

013*15 

SUM) 

OUt)U 

ft 

0 l 36 ) b 

SUMZ 

0000 

ft 

007522 

51 

noao 

R 

013556 

TB 10 

OOOO 

ft 

012575 

TITLEl 

oooo 

R 

013561 

TSTO 

Uuyu 

ft 

013637 

V Afi 

oouo 

R 

010*23 

VARCOV 

oooo 

ft 

013601 

VARERft 

anoo 

R 

013572 

WAVE 

□ 000 

ft 

01 1223 

wfbe 

oouo 

ft 

010670 

K« 

— — UOU0- 

ft— 01 3163 

WPS E 

oooo 

ft 

01 2207 

A - 

-6000- 

ft- 

0+353*- 

XBAA 

OQOB- 01360/ 

XL 

OOuD 

ft 

010506 

Aft 

□ ouu 

R 

01363s 

AHA* 

Oooo 

R 

013531 

Amin 

□o 0 0 

R 

013605 

XMX 

(1000 

R 

013610 

Xft 

OOyU 

ft 

01 2bl l 

Y 

OODU 

R 

0 l 313s 

YAWS 

onoo 

R 

013611 

Ytl 

oouo 

K 

010537 

YM 

6000 

R 

013606 

YHX 


UOuU H U13612 tT 


ou 100 

l* 

UUIOI 

2 * 

00103 

3. 

00 |UH 

4* 

0010 + 

5-* 

00106 

6* 

00106 

7* 

00 1 06 

6* 

06106 

9 • 

00 | Ob 

1 0* 

00106 

1 1* 

UU l U6 

12* 

00107 

1 3* 

001 lu 

14* 

001 II 

16* 

001 12 

16* 

00113 

17* 

uom 

Id* 

00116 

| Y* 

OOUO 

2u* 


INlT.ALUE *r> D INPUT BASIC PARAMETERS 
INTEGER Bl0ELi6l02iASEiAtigEiABE t ANsE 
DIMENSION A‘)10(6 ,6, I0l3t »£5CH*I !0, 2S I 

D I HENS 1 On SA t 10 ,2b > ,SCT I 10, 25) , XhI 25) » YH 125} ,*M2b» , aNSEI U I , 

f ttrefM 2> ,WiRC0V»2*2»,fil2),Q{2},|TtTXII2TrlT|T¥HH> 

COMMQN/P/P 

OImEsSIOn PU0D,3bl , ACUO, 16 I ,SE( 161 ,NZSE t 16) »0 En5< 16 I ,*PSE< 1*) . 

1 AU+BI Ul l6i *AE t 10> *NZ ft E t |0) , ABCS! 101 ,S IGSUO J ,AFS C >U ) , 

2 ABE t l Q ) , BG 1 |01 ,C1 r( 101 ,EScl!ui , C + I 101 , ABC I 10) t A F*M 1 q ) tABCOtlOl , 

3 «H0B2lt0).DELtlnl,SlQ0l 1 10»,B!0ELtl0l,BIOZ(ID).A«TGELUr)l, 

+ AM+fttJ-e+t 61 , *-&E l6>tAMT* J SE<6llNZA(6l,0ENABl6l,SlGB2llol UoftC t )0 » , 
b AUhSI I 61 ,S0VESC( 101 , SOVEKrI 10) , 5 CTSD( 10,25) 

DIMENSION *(2k)5),F|1l 1 ,TITlE 1112) 

DIMENSION YU01,4),YAVG(1),CV14,i),CC(4i4>,CP<4,4> 

DIMENSION 1CS< tOOl 
DIMENSION 10ATEI21 
Edl) hi AtENCE fPl t , 1 1 , ICS ( I I > 

DATA BlUtl- 7* K*,' C * i * N«,» O', •+*»,' P*,* S’l'CL*,* K * , * F E • 7 
DATA til 02 /I ,0,7,0, | 1 , 15, 16,1 7 , |Oi2t/ 

Data AMTcE L /o/.29,‘t?.7Z,i5.oJ,34.7 1 t 1 co.on,Dn.uu,au.zz,ao.uo,u0.oo, 


05120 
□ 0122 
00 122 
00124 
□ 0126 
00126 
-00130* 
00120 

00132 

00133 
0013*4 

00135 

00136 
00(37 

00137 
00137 
00137 
Q0U7 
00} 37- 
00137 
00)37 

□ 0137 
00137 
00137 
00}37 
00127 
00137 

□ 01*40 

00151 

00152 
00163- 
00)63 
00163 
UO I 6*4 
00170 
tOl 7H 
0017*4 


-216- 
2 2 * 
23* 
2*4* 
2S* 
26* 
-276- 
20 * 
27* 
30* 
31* 
32* 
•3-3*- 
3H • 
3b* 
36* 
37* 
3B* 
-37*- 
*40* 
Hi* 
H2 * 
*43* 

*4 *4 * 
95*— 


-1- 


- -OOrOo/ ' 


*46* 
*4 7* 
40* 
*49* 
50* 
5-1* 
52* 
53* 
54* 
55* 
56* 
5/*- 


OA-f A AM TRBC/. 072), .c,36B,<1671,, 2006,. 0tH2>. 0030, ,00361.00 45, .Oil*), 
i ,0031/ 

OaTA TITlEI/'CONcentrATiOn FREQUENCY 0 1 STR | BUT ION ' » 6* ' V 
OATA UTlTxtl » i !*i , I21/72M A8 SORPT 1 OM-FAC TOR * ARLA-0EN5 i TV ( 10**- 
161 / 

0 f t T A — M - T -1T t H-4-f 1-* 1-^ l 2 1 / 7 2H AC-£UMUL A TOO COOP! T — i~l~9*-*6i 

I / 

NAprEt.lST/l>lST/Cl ,Cl2,C2, AOBt AKATlO 

NAMELIST/bOATA/DEnAD.OEWSiADNbiADS.SlOS.AC, ABCS.AFS.ftPSE 
HAnetlST/OuHlO/NSTO.XM.yH^HiBG.ESC , SO VESC , SDVERR 
MAnELIST/B t iATA/A f 3C l AFB,AOCO,CR i 5ig 5l SIG02,ADNBfADB, AFBE , I COunT 

- — NAMEtrl-^f/STAT/'AH^ff, XhAX , XB ARt V A R , SDE V — — — --- — — ■ ■ 

CSCtH ■ 1*2605 

* UlST OF OPTIONS 

lPtt|«T * 0 NO 44AHEHST PRINTOUT t T VP ICAt 1 

I PR I NT * I PRINT NAhELJSTS CONTAINING INTERMEDIATE CALCULATIONS 
16 V"T e 0 WttCrttTS-ARE- INPUT- THROUGH CALtlNO A-»G«HEfiT5-l*t-/5*MPLE YARl 


ir 

1 HEIGHTS Are DEFINED AS t, in GLHFR2 

2 HEIGHTS Are CALCULATED in GLMFR 2 AS 1*/VII1 IJ/POlSSON VAR1 
* o URL V STaIISTICAL ANALYSIS of STANDARD DATA PERFORMED 

i l ype * *— i samples are red blooo cflls 
ITYPe * Z SAMPLES are biological TISSUE SECTIONS 
- hype at 3 ■■■SaHpfcreS-ftttE GEL-ATjO SECTIONS 


I WT » 
t*T * 
1 TTPt 


— c _. 


INPUT PROBLEM INFORMATION 

RE A 0 ( 5 i B ) U>aTE, HISS, man. JSET 
a F QR4+A T ) 2*5 , I5rAS,AS| 

AR 1 Te ( 6 , 20B | 1 DATE, HISS. Han » 3SET 

-CARH- l *,2A5il6»A5rAS» 

• INPUT PROGRAM PARAMETERS 
RE AO 1 5 i *4 1 (PRINT,! By 
ARIT E lfr,2 IO) IPRInT » | AT 

2 1 0 FORMAT!' LAHp 2 *,3J5I 


00174 

SB* 

L 

• 

INPUT MtCttOPRQBE PAKAMETERS 

00175 

59* 



REAO(5,7) AP,BOiCT»DT 

U02U3 

&u* 



7 format ( 4f i o,2> 

UU204 

61* 



HR! TE46,207 1 AP,BO,CT»0T 

0021 2 

62* 


207 FORMAT!* CARD 3 ’,4Fl0.2l 

002 1 2 

(.3*- 

e 

— 

— * — 

00212 

64 • 

c 

• 

INPUT sample INFORMATION 

00213 

65* 



RtADI5i4| ITYPE, N5AHP, NA tTelO 

U0221 

66* 



4 F0RHAT!3]5,FS.U> 

00222 

67* 



AH I Ttt 6, 204 ) ITYPE, NSAHP,NA,Td!0 

00230 

65* 


204 FORMAT!' CARO 4 **315,F5*0} 

00231 

69* 

■ — 

— 

T411-0— »-Tti 1 u*-TttOE^U4 

00231 

70* 

c 



00231 

71* 

c 

■ 

Input standard oata 

00232 

72* 



READ (5,6) gEL*H20,TST0,nS1d 

U024U 

73* 



* FORMAT 1 3F l(j , 5 , 1 5) 

0024 I 

74* 



AR|TE!6j2031 GE L « H 20 , T5T0 , nSTu 

00247 

75*- 

- 

*03 POR+1-AT-+*— C AfTO- 5-* t3F 1 0* 5 * 1 5 1 

00250 

76* 



TSTO * T&I o* l .0E-E14 

00251 

7 7* 



REftO 15,2) ( AE4 1 I ,061 I 4 ,ESCU » ,SDvESC( 

00262 

7 B * 



i FORMAT IA2 , 3A,3F Io,0 1 


I ) , I “ I , N A I 


-U0263-.. 

00274 
002 7 4 
QQZ7 4 

00275 
00277 

— 0 □ 27-7 — 

00301 

00302 

00303 
00}0H 
003U5 


-JSj 


«R .|.T.E<.A^au a J-4- A E 1 1 L, at, l 1 1 ,-ES t 111 , S D V E S C U 1 . I « 1 . HA > - — 
CARO 6 • ,A2,3X , 3 f 1 0 * 0 1 


aq* 2 q2 format ( 

til* c 

82* t • COHPuTE roc area density anu current intercept ratio 
03« I F t 1TYPE . NE » 1 1 Co TO 35 

•OIAGhOST JC* The TEST FOR EoUalITT BESEEM 
8-4* IT 1 !- “IrFTBtTT&lf 


85* 
86 * 
87 * 
66* 
B9* 


U UmJ U U 

00307 
00310 
0031 1 

00312 

00313 

um-I-ii 

ITU** 

91* 

92* 

93» 

94* 

95* 


00315 

96 * " 
9 7 * 


00320 

98* 


00320 

99* 

C 

00320 

100* 

C 

00320 

1 D 1 * 

c 

00320 l-tF2*- - 

— * ■ 

00320 

l U3* 

c 

00320 

104* 

c 

00321 

10&* 


00324 

106* 


00334 

107* 


00336 - 

100* 

— - 

00345 

109* 


003N5 

110* 


00350 

Hi* 


00361 

1 12* 


00362 

1 13* 


00373 

1|4* 


00374 

1 15* 


003 7 7 

t 1 6* 


U040I 

117* 


U04U5 

1 le* 


00405 

1 19* 


00405- 

1-20* — - 

— 

00406 

121* 


00407 

122* 


01)4 11 

133* 


0Q411 

124* 

c 

0041 1 

125* 

c 

004JI 

1 26*- £-- 

00412 

127* 


00414 

12a* 


00417 

129* 


00422 

13Q* 


U042H 

131* 


00425 

132* 


00426 

133* 


00430 

134* 


0Q432 

135* 



NOn-JnTeGERS May NOT he meanim&fijl. 

y>ttgo- Tcr-jrt 

cr* 1 . - EXP(“0,s/(]. 75**2) 

Ci2~n T« •* EXP(-4»S/a.75**Z) 

C2 « C i 2 - cl 

- A08 » 32. • | .01 720*C 1 * *4ft272*C2/8»)*,G0Ql/clZ/3»l4lS9 

60 TO 312 

•y rt -'Cor . -riT j uE- — — - --- 

Cl ■ l. - £XPt-4,5/( | .S*BD)»*2)1 

C12 * l. - EKPt-O.B* ( 3.7&**2i/M ,5*001 **2) 1 

C2 <■ c 1 2 - cl 

AOB * 32.*(tJ.5*C] /9* + 0 ■ 5*C2/ I 3 a 75 •> *2 - 9 . H *0*000 1 /C 1 2/3 . I 4 1 59 
312 CONTINUE 

— a-r at i u » -i¥tiyetg - 

(HI I TE ( 6 1 D 1 ST ) 

35 CONTtNUE 


this Loop forms atomic 

'•P & frS-f T-y AftO-A a-StHtP T I ON 

each standard 


NUMBER Aro SYMBOL ARRAYS, ANO COMPUTES fl REA 
COEFFICIENTS FOR EACH- AfTA-tYZEO ELEMENT }N - 


DO HI NO * l,NSTU 

REaO|5«S| NAO i I A st I H ■ AMTASE ( I I , I » 1 , r J A D l 
5 FORMATtlt ,3 a,6IA2, 1X,F6.0| 1X1) 

«1TTi 1 ?t-6-r2li5T*+TAD-, t SSE t 1 1 i AttTASE I r 1 • 1* 1 ,NADl 

2nS FORMATS' CARO 7 ■ t I l , 3 X , 6 < A 2 » 1 X , F* . 0 , t X )) 

IF (N3TD«E0. 1 I GO to 321 

RE AD (b, 305 1 IPS. IANSE t II tSCTI 1 ,ND) i 5CTSD • I ,NDl , I«1 ,NA» 

305 FORMAT I l2.2K|RtA2,lx,F6.0>lX.F6.Oi!XII 

A«lTEt6.2U6J IP5,tANSEtI>,SCTI|, N 01,SCTSDU,NDIiI a, l*NAl 

"208 F O R MA T t -*-~CTt Rp- S ♦ f J 2 , 2X , 4 t A2 . 1 X , f6 . 0 , 1 1 1 F 6 . Ol t* H 

DO 32 I * | , N A 

IF t ArjsEl I 1 .E0« AEu » 1 60 To 32 

i>RI TE 16 ,241 ] NO.NSTD 

2ql FORMAT! 1H0,5X| 'ANALYZED Element CHECK ARRAYS AN5EIII ANO AEJll, I 
1 * 1 ,nA. ARE NOT IDENTICAL. 'i/.SX, 'CHECK CAROS 7 AND 8 FOR STAnDAR 

GO TO 999 
j2 CONTINUE 
321 CONTINUE 



■ construct a i dm i c number arrays nzaiij 

-*WU~*H- A L Y2 EtJ-Eb£HE N 1 S 
IF 1.411. NL, 1 | GO TO I l 2 
UO l I I * | i N AO 
00 12 J * I i 100 
IFIlCbtUj.NE.ASE(l) J GO To 12 
NZAtlt * J 

-e g - T o Lt 

1 2 CONTINUE 

IF1N2AI I J *NE. 01 GO TO 1 1 
n«ITE ( 6 , |U l ) I 


ANO N2AE1I1 FOR ADDED ELEMENTS 


QOH36 
0QR37 
00111 
UQ111 
Q011 7 
UD1&1 — 
UQ1S2 
UG1b3 
0015S 
0015 7 
00162 


1-34*-. 


-Vb l-F-OR ^ - T -t-t- K lJ , i-QJ-r-^titOEO Efc^HEWT SYMBOL » r I VX*t*AS--N07 -EOUUd ! ) 
137* GO To 999 

1 3&* li CONTINUE 

*35* 00 13 1 * 1 i HA 

1 HO* DO H J «i,|QO 

111* IF( lc5tJl.NE.AEt 1 I) GO TO (1 

-+»** — n^ao-i-h— *- o - - — ■ 


113* 
111 * 
»15* 
ll6* 
l *i 7 • 


® T ° 1 

ti COuTiuu 


3 

UE 


IF (NZAE I J I .NEtO) GO TO 1 3 

*R-ifri*nuzi i 


OB-163— 

00161 

- M-a* — 

1 19* 

13 

-QO To ftT- 
CONT I NUE 


001*1 

1 50* 

C 


U0161 

151* 

C • FU»iM SYMHOLi 

OENSI 

004*6 

152* 

112 

CONTINUE 


00167 

153* 


DU 2| I > 

l » 10 
rEL* ( h» 

001-72— 

I’S 



00173 

U0171 

00175 

00-476 

00500 

00*03 

U0SU1 

UQ&Ub 

00507 

005)0 

00513 

0051* 

005*11 

005*2 

005*3 

00521 

005*5 

"005*6 

005*7 

00530 

00531 
00531 
00531 
00533 
00531 
00535 

00510 

00511 

00512 
00511 
0C5N7 
00550 
00552 
00552 
00552 
o05*H 
u Jbb7 
00560 
oObb 1 


1 55 • 
156* 
157* 
158* 
159* 
"1-6TJ*- 


1 6 1 • 

I 62* 

I 63* 
161 * 

1 65 * 
T-fr6* - 
) 6 7 * 

1 66* 
l 6 9 * 
l 7 tl* 
17 1* 
l-wt»- 

173* 

171* 

175* 

1 76* 
177* 
178*- 
17V* 

1 au* 
let* 
l a*» 

1 S3* 
lot* - 
1 B5» 

1 1)6* 

1 B 7 • 
fee* 

1 S9 • 

1 vu* 
1V1* 

I 92* 
193* 


5£ ( [ t « BlOELI I I 
NZSEt 1 ) * 8102(1) 

DEnS(J) * DRYGEL*AHTGEL( 1 ) * l *0E-02 
21 CONTINUE 

DO 22 1 * | ,NAO 

f t T ’J i/rt-f) 


OENAjin » AHTaSE I I) *P(N I U • l .0E-06 

22 Con T I NOE 
K » 10 

00 23 1*1 ■ MAD 
00 24 J * liK 

— ■ E^t 5 t£+-J n G 0 TO 25 _ — 

2 H CONTINUE 
K « K+l 

NZ5E (K> » NZ A 1 I ) 

5E ( K ) '* AbE(I ) 

DENS f K 1 - DENAD t I i 

— _ 60 — f-fT~2-3 — — - . 

25 CONTINUE 

DEnS ( J) * OENS(J) + DENAD t | ) 

23 CONTINUE 

c 

c « calculate akea density OF Dk y standard 

irinr~*~-K 

AOS * 0* 

00 5 1 J * ( ,NSE 
AOnS(i) • DENS ( I ) *TSTU 
AOS m AOs * ADN5 ( I ) 

51 CONTINUE 

rtPSE(I) * 100 «*AdnS 11)7ADS 

52 CONTINUE 

IF (NO * NE ■ 1 ] GO TO 113 
C 

C • COMPUTE ELEMENTAL mass ABSORPTION COEFFICIENTS FOR analyzed lines 

DU If I * t iNA 

N * NZAEtl ) 

AA * P t N , 6 1 
00 H I J « l ,NSE 


L$-V 


00565 
00570 
□ 0571 
00570 
00575 
— 005-7-6 
00577 
00600 
00601 
0Q6Q4 
00605 

- 0060 

00607 

006)0 

00613 

00614 
006)5 

—006-)?- 



00 43 M - 9,17 
SDG&* Pj NZ ,H) 

If (VlAVE.LT.EDU£> 

43 CONTINUE 

COh ■ P<N-d,3D* 

— EX » 2.22 — ■- 

GO TO 45 

44 CON n PlN2,rt*12l 
IFtN-10146, 47,48 

46 EX- ■— P ( N 2 I 1 9 I 
GO TO 45 

— q7 Ex ■ f Ut z r 2 tr) 

GO TO 45 

4« TPt«»l3l 4? ,49,50 
4? EX * 2,60 
GO TO 45 

50 IF [M.LT. 1 7 1 GO To 
• EX- 


1.0 TO 49 


49 




00620 

2)3* 


45 ACt|,J) « CON • 6 A VE * * E X 

00621 

2)4* 


41 CONTINUE 

00621 

2)5* 

C 


0062) 

2)6* 

c 

• COHPoT-E MASS ABSORPTION COEFFICIENT 

30624 

2)7* 


1)3 CONTINUE 

00633 

? i 



00630 

219* 


U U — 9‘J ■ 1- — »- f ftHr 

ABCS ( J 1 * Q, 

00631 

2 20* 


01) S3 J • 1 ,NSE 

00634 

221* 


PABC ■ AC l I , J I «AD«5 ( J 1 

UQ635 

222* 


ABC S | | )- * ABCS • I 1 * PABC 

00636 

223* 


53 CONTINUE 

— 00636 — 

— 224-*- 


. 

00636 

225* 

c 

* COHPyTE ABSORPTION FACTORS FOR ANAL' 

00641 

226* 


00 54 1 « 1 ,NA 

00644 

227* 


POE * A0 CSU)*CSCTh 

00645 

2 20* 


AFSUl * 111 - ExP 1 “POE I 1 7P0E 

UQ646 

229* 


54 CONTINUE 

-00646— 



— f 


00646 

23)* 

vr 

C 

» COMPUTE area DENSITIES for analyzed 

0Q650 

232* 


00 Ss 1 * 1 ,NA 

00653 

233* 


DO 5* J a 1 ,NSE 

00656 

234* 


lFlNZAEl 1 I.EOiNZSfcTj) ) GO TO 57 

00660 

235* 


56 continue 

00663 

236*' 


st-c-ontinue - 

00663 

237* 


SIGS(I) - A0N5 1 J 1 

00664 

230* 


Sb CONTINUE 

00666 

2J9* 


00 1 1 H 1 * 1,NA 

00671 

240* 


SMItnOI * AF 5U) » 5 1 G5 1 1 ] 

00672 

24 I • 


l l 4 CONTINUE 

00*74 

242* 


- I’PtFPR I Nj 1 NE-f O-}— rtrlTE 1 6 i SUAT A | 

007UO 

243* 


IF (NSTD.EN. 1) GO TO 116 

00700 

244* 

c 

■ LOAD ANALYSIS INPUT-OUTPUT MATRIX 

00702 

24b* 


00 3 | t * l , N A 

00705 

246* 


ftttiOt t 1 1 »Ntrl » AhIASEIIJ 

00706 

247* 


ANIOI I ,2, Hoi » I0U,*S|G5( U/Aos 

00707 

240* 


AN ttm-T3y+ttrr *— SttiS ( T 1 

U0710 

249* 


ANIOI | ,4 ,ND I “ AfS| | 1 

0071 l 

2bU» 


31 CONTINUE 

00M3 

25 | • 


IF ( tPS.NE i 1 ) GO TU 111 


ENT OF STANDARD FOR ANALYZED LINES 


LINES IN STANDARD 


ELEMENTS IN STANDARD 



0U-71&-- 
007 1 6 
00716 
00716 
00716 
00720 
-U0733-- 
UD726 
0072 7 
00730 
007 32 
00734 
-00735- - 
00736 
007*1} 
00742 
007*13 
00744 
-067-46— 


-25 2+ • 

253* 
254* 
255* 
256* 
25 7* 


- iND*" 


— N [y- 


l l 1 CONTINUE 

• COMPUTE estimates OF STATi 
AND COMPUTE MOgT PrOqaBLE 

DU US I ■ 1,NA 
H fl-1 ' 17 N p * }iN *Ttr- 


5 T 1 C A L PARAMETERS FOR 5ET OF STAnuAR^S 

count for the principal standard 


00747 

00750 

00751 

00752 
00756 

-00766- 


259 • XrtlNDJ ■ SX [ I i NO ] 

-260* - YMUNOt *-5 cT|I«Hd) 

26 | * IF t I «T*Eq*0 > tYMIl) ■ 1 «D/SCT50C 1 ,ND) **2 

262* 11/ eOHT-TwOE— 

263* HD(M b 25 

-2*-H-*— — C A tTfc- OLI If^^lfrTDTXM. rM-FVfH-T-HD lH-i IttT r£OEF rV-*fteftR-rY-frM5oV3- 

265 6 DO 34 ND*l*N$TD 

2*6* £5C«X(lrN0J * COEFtt) + C0EF121 SXII.ND) 

267* ANJOf 1 ,5, NO) * Scfll.NP) 

2*8* AN TOT l-f 6 , NO 1 ” EsCMXtltNO) 

267* 34 CONTINUE 

27 0* Bu r 11 * -CNEf-M-F - - - •*•• 

ESC I I } » COEFUI + C0EFI2 » *SX I I , IrjPX ) 


271* 

272* 

273. 

274* 

27b* 

-2-76*- 


50VESCH1 = SURTUBSIVARCOVUil) 
SPvCRHU) * SURTtrtBStVARERRI > 

iPriPRiNT.NE.ai ivr i te i a t sen i o i 
punch 2, i aei I i »bou j ,esc( i i »sdvesc ( l n 


VARC0Vl2,2)*5X*I,lhDX |t*2|l 


00756 

277* C 

• MICROF ILH ^TauDARO 

COUNT 

007*4 

276* 

XMX * 0. 


00765 

27 9* 

TUX * 0. 


00766 

280* 

00 lie J » 1 iNSTo 

0077 l 

2B| * 

XMtJJ * X M i 0 I • 1 ■ 

Qt + 06 

OQ772- 

-26-2* 

— Ytt-t-P I *- TM-t Jt-*t • 

OE-03 

00773 

203* 

1 F < XM l J) .GT. XMX ) 

XMX - 

00775 

204* 

IF ( Y M t J I ,5T *YHX) 

YMX » 

00777 

2 (ft,* 

11 0 C OUT I NOE 


01 0U| 


XL ■ Q. 


01002 

£(./* 

XH ■ xMX 


01 063- 

2e»-« 

Ttj“*-tn 


0)004 

209* 

IF iCqEFI 1 1 ,LT»0. 

I YB * 

01006 

290* 

YT * YHX 


01007 

291 » 

Qi 1 } * 0. 


010)0 

292* 

ai-2) * XhA 


uian 

293 • 

Rill - COEF l 1 1 • l 

.UE-Q3 


vS. ABSORPTION-FACTOR • ARE A— DENS | TT 


0-1012- 
01013 
0|014 
0)015 
01020 
01020 
01021 
U 1 023 
01025 

01025 

01026 
0103 1 

01034 

01035 
O}037 
01037 



XNUI 

YHIJl 


* COEF ( 1 ) • 1 . 0E-03 


2F*lw0*-09-*XMX — 

CALL 4UIKMP[-:,XL»XRiTB,TT,lH*,lT|Tx,lTITy,NSTD)XN|YM) 

CALt OOIKMLIOiXL.XR^OiYTiIH iITITX,ITITY,-2,Q,R1 
6 RI TeI 1 / , 104 | NZaUI 

[04 FORtlATl 1 h , *-*3BX I * ACCUMULATED-COUNT vs corrected AREA-DENSITY for el 
1EHENT 61th 2 * '.12) 

Ib-COTS-T-Itjoe- — - 

WKITei 6,2U) 

2 1 t FORMAT 1 1 h 1 » * 

I standards * | 
oo a VI.*.*- *.*; 

MUTE (6.212) 

?l* -F t > 4 M A T t l -M O i ' P/~ri* ,*2 ( i 
AftlTEl6,213f 

2 J 3 FORMAT ( 1 MU . • NUMBER * i 6 X t • CONCE n TR aT | ON * , 2 X , ’ 


STATISTICAL ANALYSTS Mt) C^O^ROBE 

.Ci» •• ” ^ 

. ^ ' V A-U» “ 

A£ I I J 

15 THE EtEMEHT ANAtVEEDM 


*«? - - 

OATA 


FROM A SET OF 


WEIGHT 


l » AREA 'jaXt' ABSORPTION* |5X > »MFASuRE0 , t7X 1 *ESTlHATLD , i/jlX. 


VioitN-i > 


~*ir 


M-0^7- 3 (Q » 




LI (037 
0(040 
0(043 
0(043 
0(054 

: — G405-4- 

OiOSb 

0(057 

0(065 

01065 

0(065 

0 (-0-6-5— 

0(056 
0(070 
0 1 07 ( 
O(07( 
0(07( 

0(071 

0(07 ( 
0(071 
0(071 
0(073 
0((0( 

U-H-02 

0((05 
0( l 07 
011(3 
OKU 
f Q(U3 

t-i Ol-t-i-4 

10 0(116 
oi ; 21 

0(127 
Q( 1 27 
0(127 
-0-1+27- 


3(1* 

J-12*- 

3(3* 

3(4-*- 

315* 

3(6 * 

3(7* 

3 1 66 
3(9* 
320* 
32 ( * 
-3 22* 


FACTOR KBX,* COUNT 


f-*-t*-x-, 


3 • FACTOR KSX,* COUNT «i7X,i COUNT ') 

00-9-4 ND = ( f NS TO 

fJR]TE£6,215) ND • AN| 0 1 1 i l »N0 1 i AN| 0 t I »2 »N0) * AN(0 ( ( ,3 ,nd) , 

— l-Atr(0-M-T^ jfJ03rANl0l I ,5 in 0) , ANln< ( (6 lNO» 

2(5 7ORl1AT(lH0i3X,13,9X,3X»F5»l*7X,tX(E8.3,6X,2x.Ee.3i5X(3XiF6*6(6X» 


323* 
■324* 
325* 
326 * ■ 
327* 


94 CONTINUE 

6RiTt{fr r 2l61 B G I I 1 , CT.SDVESCI () ,n 5TD 
2(6 FORMAT (1H0, l X<F6,0, • (S THE ESTIMATED 
i F-4-»£)-(*- -SECOND COUNT 
2ANDARD DEVIATION (OR 

09 CONTINUE 
1 T 6 COwTrnUE 

IF (NSAHP.EO.D) 


RESIDUAL BACKGROUND FOR A 
TlHE.»,//,lI(,F6.D^ , IS THE- ESTIMATED ERROR ST 
INTENSITY COUNT ON* ,/>BX( 'THE PRINCIPAL STAND 
I*T- — ~ 


GO TO 999 


329* 
330* 
331* 
332* 
333* 
-3 3 - 1 4* — 
335* 
336* 
337* 
338* 
33 9* 
-34 0 * 


0)130 
0(131 
01 (34 
0 | (35 
01(40 
0-H-41-- 
01(42 
01(44 
0(147 
Ol 150 
U ( 1 52 
0| (52 
0(152 

01153 

01154 
OK 54 
0(155 
0((56 
011*1 
U( (62 
01 162 


34 I* 
342* 
34 3* 
344* 
345* 
-346-*- 


347* 
34 0* 
3 4 9 • 
350* 
35|* 
-3&2*- 
353* 
354* 
355* 
356* 
357* 

- 35-8 
359* 
360* 
361 • 
362* 
363* 
-36-4*- 
365* 
366* 
367* 


-*— F&ttPuftiT-A-NArfcrS-r^— Or ONE-SftHPLE EACH T-TME 'THROU G H L OOP- ■ -- 

• CHECK SYMBOLS AND STMBOL ORDER OF COUNT DATA 

REA0L5|3) ( ABE ( ! K | K 1 >NA 1 

3 FORMAT 16 l 0X ( A21 ) 

-DO— fr| — i— = — TlNA — • — - — ■ -— — 

IF IAE11 ) .Eg.ABEl I 1) 50 TO 61 
KRtTEL 6,1031 NA 

J03 FORMAT I IhOi 1 OX t "CHEMICAL STMBOL OR SYMBOL ORDER IS IN ERROR FOR N 
+A « » ,12) 

GO TO 999 

— 64 -- eOfrT-I-tttTe -- - — 

DO 99 NS = KNSAmP 
■ REaD(5t581 { CTB 1 1 1 t | * 1 i NA 1 
56 FORMAT (6F 10.01 

* form Initial density array for biological elements, 


nbe ji hse 
(FtltYPE - 2(63(64,65 
63 CONTINUE 

DO 66- ( * 1,(0 

ftFaEui - amtrscki 


66 


69 


CONTINUE 

DO 69 I -■* l,UA 

CR i I 1 * CTBt 1 1 / t ESC ( I I 

CONTINUE 

GO TO 67 


BG (I ) > 


• SPACE LEFT 
64 CONTINUE 
GO TO 67 


HERE FOr COMPUTING TISSUE DENSITY ARRAY 


65 COnT I HUE 

Do— 62 I * I i 'N A- 

CR l 1 1 * (Ctu [ 11 ■ 
62 CONTINUE 


•bgi i ) i/ttsci n-BGt l ) ) 


c 


-40 


0 mi- 
ni 167 

01170 
01172 

01171 

01 171 
1*1173- 
01171 

01175 

01176 
0120 1 
01202 

0 1205-— 

01 206 
01207 
01207 
01207 
01212 
UJ21S- 
01216 
01216 
012 16 
01217 
U 1 2 1 7 
01217 — 
01220 
01221 
01223 
0|223 
01223 
01221 
012 27 
01232 
01231 

01236 

01237 
01210 
01212 
01211 
U I 2 15 

01250 

01251 
uta&t 

01251 

01253 

□ 1255 
u 1260 
01261 
0 1 263- 
01266 
01270 

01272 

01273 
012/1 
01277 
U1 JUO 

□ 1 3U2 
01303 


3itO+ - 
3&9» 

3/0* 
37|* 
372c 
3 7 J • 
-3-71*— - 
375. 
3?s* 
377* 

37 B* 
379* 

3 bo* ~ 
3a j • 
362* 
3a3* 
3ai* 
305* 
-3Bt5t * 
337* 
336* 
339* 
390* 
391* 
- 39 - 2 *- - 

393* 
391* 
395* 
39ft* 
397* 
3 yo*- 
399* 
100 * 
1U1* 
102 * 
103* 
101»- 
105 • 
106* 

1 07 • 

hob* 

109* 

110* 
111* 
112 * 
113* 
111* 
115* 
1 l 6 * 

1 1 7 * 
118 * 
111* 
120 * 

12 1 • 
122* 
123 * 
121* 
125* 


-0 0 '6 ( 1— 1 - 


— ItNOE 

AOmBiJI * A0N51 IJ-TBI0/T5TD 
60 CONT 1 ItUt 

ADB * ADS*T310/T5T0 
67 CONTINUE 


-T-tO0PtrTe"H**5~(rB5OfipM 0N COEFF* OF 310, S AHPfcE F&R- tl Nes— A4i A t Y Zfc D- - 
ICOUNT = 0 ' w 

7Q continue 

do 7 1 i a t ,NA 
ABC l|) * O, 

” DO 7 1 J * t iftaE 

— -PTSB^"* r ^TCt-fnri-*-AoiYa t j> ■■ — — - . 

ABC! I t * A BC f I ) * P A B C 
71 CONTINUE 

* COMPUTE ABSOHPT 1 ON MCT0K5 FOR ANALYZED LINES IN BIOLOGICAL SamPLE 

DO 72 I s 1 ,nA 

■ POE—'* A BC 1 1 ) *C5 C Trr • .. .. 

AF B I 1 ) * 1|« * ExP ( -POE ) ) 7P0E 

* COMPUTE TH E ABSORPTION CORRECTION ARRAY FOR THE ANALY^eO ELEMENTS 

AEJCOlr ) = AFSI I I / AFB t 1 ) 


-e 

c 


¥ COMPUTE AftTP-OttfS-tTY' ARRAY OF ANALYZED ELEMENTS -f-ftf" BMHrOfrtCAL SAMPLE 
5 | 0B2 j I ) * SI6S1 I 1 *C« l 11 * ARAT 1 d*ABCO III 

72 continue 

KoUaiT n ICOUNT + [ 

• UPDATE area density AND AEtGHT FRACTION ARRAV5 OF BIOLOGICAL SaMPLE 

” 170“ ?S~ I-e— I ,WA 

00 76 J a 1 , nBE 

IF (N2AE1 1 ) .EO.NZSE [ JJ 1 GO TO 77 
76 continue 
7? continue 

AONOI J) a S1GB2I I I 

•••ys - e offT T nu e— — - .. 

iF(lTYPE.tO, 3 ) A 01 *B I 8 )* t A 0 NB 1 b) / 72.997 + ADN 8 t 9 I / 3 H 1 UO) * 35 . 1&7 

ADB * Oi 

DO 7 8 1 a J ,fJBE 

ADO a Al>B ♦ AONB ( 1 ) 

7B CONTINUE 

• test for Reiteration criterion 

IFUCOUNT'LE .2) GO TO 73 
Oo 71 1 a l ,NA 

deli I ) * <S IGB2I I 1-SIG01 < I I ) 7SIGB2I U 

71 CONTINUE 

DO — SJ— T — sr J jNA -■ - . .... ... 

IF IAbS1Dee{] 1 1 ■ G T • 0 ( 'J 0 1 •And*1C01InT*LT«101 GO TO 73 
33 CUhT I HUE 
GO TO B1 

73 continue 

00 00 I a | , N A 

^ tg u - n it-* 3 lira 2 1 tt 

au continue 
GO TO 70 
01 CONTINUE 


-41 


£- 


ajiUM * 126 * -t-R+HMit T Htre4-Go-T0 »6 - - 

U|JU6 427* OG B 2 I b l,h)A 

U| 3 1 J 428* N b M2 A ( I I 

01311 42 9 * t 

0I3U 4J0* C • COMPOTE THf YOLUpE DENSITY OF ANALYZED ELEMENTS IN 5ELATIN SAMPLES 

01312 43 J * RHOBZ(I) = S I GB2 ( 1 ] /T0 I 0 

01-31*3— — 4-32* CONC 1 1 J — »— (+ttO*2t l ) + t rUE+08/P <N i n- — — 

01314 433* 02 CONTINUE 

01316 434* i )6 CONTINUE 

01316 43s* C 

01316 436* C * F URtV; if£ I GHT PERCENT MATRIX 6 TPM A ( N s AMP » N A 1 

0131? 437* DO 7? I c 1,N0E 

0-1322- 43S*- tS-F-p-E-t- M - -a&NO-t 1 1 / ABB 

01323 439* 79 CONTINUE 

01323 440* C • LOAD ANALYSIS INPUT-OUTPUT MATRIX 

0 1325 44|* DO 88 I = l ,NA 

0133U 442* ANiO(Iil,NS( «* CTtJ(l) 

01331 443* ANlO<Ii2,Msl = A p u { | ) 

tit-332- - 4N4»- — — — AN-t-0-r[-r3TNsf— *t- -ftcto t I I - 

01333 4 4S* AN|0(l,4,liSl tt SluBZIl) 

01334 446* ANtOMibjNSl - I (JO , « S f GB 2 U ) ? ADB 

01335 447* A N 1 0 ( | * 6 , u S ) » CO«C I 1 1 

01336 440* 00 CON-tMJUE 

01340 449* lFtNS,5T.2l GO To 9* 

O-ta** 4*e* 1-FrtPTTHirTTJE-FO l- VNl TE I6,00ATA) - 

0|346 45|* 99 CONTINUE 

0(360 4S2* I P | 1 TTPE .to . I > GO TO 40 

01352 453* NKI T£ 16,22 l 1 

01354 454* 221 FORMA T I I H 1 , • ANA L YS|5 OF &ELA T | N-B A5E SAMPLES’! 

01355 455* DO 87 I * 1 ,NA 

01360 456* - NR - [ Tg 1 6 x 2-22 > -ABE ( 1 1 ,B 6 < 1 I — 

0*364 457* 222 FORMATUHU ■ // »30x , A2 , ’ ]S THE ELEMENT ANALYZED' > 

01364 458* 1 / , 30 X , 1 5 , ’ IS Tut BACKGROUND COUNT'] 

01365 459* AR1TE16,223] 

01367 460* 223 FORHsT I lhl»i ' SAMPLE ' ,BX,' INTENSITY’, 6 X.' ABSORPTION • i 5X,'ABS0nPTI0N 

U1367 461* area *,Sx,i AE I GHT ' « 5 X t ’CONCENTRATION’ */, 1X,* NUMBER', 

01367 462* Tf-s-x-r’ — eourrr — *r?x i * factor * 1 5x*‘ correct-! on’iSx* ‘density’ *sx i 

U 1 367 463* 3 ’PERCENTAGE’ *5X , « IN hEU/L »J 

01370 464* DO 07 NS * |,ftSAKP 

01373 465* AHlTEl6»Z24! NS , A N I tl < I , l > rt S > i AN | 0 I ! , 2 , NS 1 , AN I 0 I I , 3 , N S > i 

01373 466* 1 ANTOU .H»N5I ,ANiO| | ,5,NSI ,AN10t I ,6,H51 

01404 467* 224 FOHHATIIh , I X , 2 X , I 3 , 1 0 X , I X , f 6 . fl , (j X , 2 X , F 6 1 5 , 7 A , 2X , F 6 , 5 , 7X , E 1 0 . 6 , 5 X , 

31-404 468* - 1 - 2x , F 4T-2~j-Vy ppS-« 1 1 ■ 

01405 469* 07 CONTINUE 

□1410 470* GO To ?99 

01411 471* 9* CONTINUE 

01412 472* ‘R 1 TE 1 6 , 23 1 S 

01414 473* HO K = 1 , N A 

01417 4/4*-- - 20’ , N*+6H tTfl * lfil /» ,4SX ,‘8 UAnT1 T AT! VE ANALYSIS OP* RB 4 ' SE?*-’ 5< / ) , 4QX , * 

01417 47b* .. INPUt-UUTMUT I ABLE » , 3 I / ) . 40X , » 2 , CONCENTRATION FREQUENCY DIST 

01417 476* *HI(5UT ION* i3 ( /» »40A t • 3. CONCENTRATION CUMULATIVE 0 I 5 TR I BUT 1 ON • , 3 I / 

MH17 4?7* • I , 40X , ' 4 . CONCENTRATION CORRELATION ANALYSIS’! 

01420 470* fttTI T E I 6 i 2 3 Z I ABE(M 

0l\‘23 479* 233 FORMAT I lH 1 , ///,30x ,a2 , * IS THE ELEMENT ANALYZED*) 

J 1 424 480* - VHt-l-Tet 6,233 Y - 

U142t 431* 2 33 FORMA T t 1HU , ’ S A Hp L£ ' « 8 X,* INTENSITY’ , 6 X. 'ABSORPTION' ,5X**. ABSORPTION 

01426 402 * I ' , 8 X , ’ AREA *,5x,» AET5HT 'iS*i/*]X,’ NUMBER*, 

01426 403 * 2 8 X , • COUNl 'i7»,i FACTOR » , S V , ’ CoRrEC T 1 ON • , 8 S 1 ' DENS I T Y * , & X , 



4-42 


U ^ ■ — H ■ — 

01427 4B5* 

U | 442 4-B6* 

U i 4 42 4 6 7* 

0(442 4S8* 

0(442 4B9* C 

‘ 3)442 - 490 * 1 

0(442 491* C 

0(444 492* 

0(447 493* 

0(450 494* 

U(4b2 495* 

liHil -496-* 

0)464 497* 


0 ( 9 b 7 
Q(46( 
0)463 
U ( 465 
0 (-466 
0(467 
0(472 
0(473 
0(474 
0(476 
OH77 
0 1 500 
OlbUJ 
0 (502 
0(503 

0 l 5U 4 
01507 
oibi i 
0(511 
0(511 
0(511 
015)3 
015)6 
0)52) 
0(522 
0(525 
0(526 
0(530 
0(53( 

01533 
0(533 
U ( 533 

01534 
0(537 
0)540 

01543 

01544 

01546 

01547 
01550 
0 155 1 

01552 

01553 

■J 1 5o2 


49a* 

499* 

500* 

50 i * 

602* 

503* 

504* 

505* 

506* 

507* 

-boa*-— — 
509* 

5(0* 
b l 1 • 

512* 

5(3* 

- 5 -14*- 
5 l 5 • 

5(6* C 
b) /• C 
516* c 
5(9* 

52u* - 

52 I * 

522* 

523* 

524* 

525* 

'526'* 

527* 

528* 

529 • 

530* 

53 l * 

532*- 

533* 

534* 

535* 

636* 

537* 

63a* 

539* 

640* 

54 1* 


- .... — 

00 97 NS * 1 i NS A hP 

**+J TE16 .234 1 N5 1 AN I 01 K t I > NS I . ANj o ( K , 2 . NS ] , AH 1 0 IK , 3 ,N5 I » ANi 01 K , 4 .HS 
») tA"JI0IK,5,USI 
47 CON T JtlOt 

•r-coitftt-Tr-5-f7fnw-*r parameters and pr-jnt cumul a - t - i v e- d istr t su t id h aho 

HlS’TOGRAM OF SAMPLE SET DISTRIBUTION 
00 9( 3 = 1 i NS A HP 
X ( J ) * AN 10 IK *5 (J) 

91 CONTINUE 
XttA* * X t 1 ) 

* nl r r »tu ■ ■■ — — 

00 92 I = 2.N5AKP 

IF 1 A { | > ,G1 . XHAXI AMAX * X t I » 

(FIX ( l ) tLl . XHIN1 XHJ N = X 1 I > 

92 CONTINUE 
sum! * 0, 

SUM 2—*- O* — — 

DO 93 I * J , HSAHp 
StlMl * 50H l + X I | I 
SUH2 » SUH2 + X I | ) « X ( 1 I 

93 CON T ( NOt 

X8 aR = SuHl/HSftHp 

V rtrt -g — (-5 UH-2 HSaHP * X 8 A P*X9 AH 1 / 1 t;5 AHP-* 1 I ----- 

SUE V » SQRTIVAK) 

Nfi •* NSAHP/5 

C«LL H 1 5 ] tX,N5AMP,NG,XH|N.XMAAiF,0.TlTLE|) 

CALL 161VElX t NSAMP,N6»XHlN,XMAX,F,| 1 0l 
ARlTEli.STAT) 

-96-eUlr-M-MOt— - - ■ 

IF I NA,ttl. 1)00 TO 999 

* PF.RFOR* CORK£LAT 1 0 1| ANALYSIS BETWEEN ANALYZED CHEMICAL ELEMENTS 

DO 36 I * | i l)A 

Y(Kit) “ An 10 11 ■ 5 i K ) 

36 CunT I HUE 

CALL COftANlYtNSAnH|NA ( Q|UiYAV(i,CViPCORtl01t4J 
*R1 TEl6t225) 

225 FORMAT I 1H1 , ///,3ox, *CON c£NTRAT ION CORRELATION AN AL Yj i 5 * > // 7 ) 

— Hfft i T EtfrYagyr — — 

226 FORMAT l t HU , lux , * CHEP ICAL 1 t I ox , « VARIANCE * i I OX , • CORBEL A T t ON * , 1 (1 X , 

1 ’CORRELa I f ON *,.' 1 1|X, , ELEHENT5 ',iQX,*C 0VAR1ANCE*,1DX, ’COEFFICIENT* 

2 , iox i 'probability * j 
00 3 7 I * 1 I N A 

J A = 1 

-fitt er it-* jAyNA . . .. 

CCt I *= CV 1 1 I Jl/SIJKT lev 1 1 . 1 : “CV I J, J M 
Ml I «EQ. J) 50 TO 3 0 

CP t I , J ) * 1,0 - RCOHftEICC I I «•) 1 (NSAMP 1 
GO TO 39 
36 CONTINUE 
- CMth JT * 1*0 
J 9 COnTInUE 

WRITE (6.227) AflE [ I ) , ABE ( J I , C V l I , j ) , CC ( | , J J , CP II . J 1 

227 FOHHATl iHt), >0X,A2,t - « , A 2 , F I 9 . 9 , 2F 2 t . 4 ) 



a m* 
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iH5« ... 

0iS66 bH3* 999 STOP 

OU »*7 S»*H« EH£J 

E.NU Op COuPILATlON* 1 y [ jujNOST I CS • 

BJcHAP SYMBOLIC 

~ ~a*Ott*p — eotrE ffetoOrf*ett - ~ 


17 NOV 71 lc:20;22 o 

1^-fltrv— 7-1 Htt 504 22 - — J-- 

0 


0 ISB03H2 

lH 

5*H 


*1 

t 

0 1 567242 

U 

230 


1 deleteo ) 

lOEcrTCflJ 



- ui ,F0R + *.-J14iUlG)C-t-U f vt04i. 

UN 1 V AC l 1 L) a FOKTRAN V LEVEL 22U6 QU29 A »<EXEC8 LEVEL E120I0009A} 

THIS CoHPJeAT [ON«aS OUnE ON IV nOV 7! AT I0»2i:03 

block data 

sturage used: code.ui qqdooui dataiuj cocoon blank cohhon[2i odoooo 
cdhhor Blocks! 

0U(j3 P Q07020 

STORAGE ASSIGNMENT (BLOCK-r TVPE» RELATIVE LOCATION, HA HE I 
OlluU I QOOOOtr J * 00l>3 ft OQOtlOO P 


block data 
etmrtON/P/P! lOOl 3*i 

Data [ p t 1 , ji , v»i , 36)/ i,oos,2hh, 
-* T&0QT rMOr rOeO(-95 L 9«OeO|» , 79 

* ,000, ,000, <000 , •DOG, 

♦ , DO , pQQi *000 i <00 » 

• <00 ■ <qq> •> oo i <uo, <on, 

» *T0G9033tE01 , -.5 108B296E-2 , 


■ < | | 88669SE“3 > 


SB7B323SE-5, - 

)6>/ ^<003 , 2HWE , 

* <000» <0Q0, .000, 999*000,999 

* fOtjOl ,OUO, ,000, *00D| 

* , ilU • <00i < OOq , <00, 

► too I ,00 1 < 00 , *UQ| * 00 , 

* • J00&97 lOEOl , «.72iaA7QIE-2 , 

-» ■ uill|m E *3i — , Aft»fB2?0£~5 , - 

0 A T A I P l 3,3) , J" 1 ,361/ *,939,2 HlI, 

*22b *G[)0 , »00L , ,000,226*600,999 

* <000, .000, .000, *000, 

»2*&B, 2* 7j, *135, *00, 

• <00, >00, <00, « 00 r *00 , 

rT&0*S^irT£t>tT "it 026 ?098E“t * 

• »»19H9a32lE-3, .9195OO60E-5, 

D A T A ( P ( 9,J1 iU"t ,36)/ ?»0 1 2 , 2HRE , 

*tis*noo, . 000 , .ooo,ii ><oon, 99 ? 

• *noo, .Qoo, .ooo, *ooo, 

■2.06,2*73, *350, *03, 

-*■ ,-Q-CK rOO-r- ,00, • DO i »00» 

• • 1 ill l 8V9QE01 , - , 1 H2S5909E-1 , 

• »«2Bi 15U22E-3, , 1 J93761 1E-H , - 

DA T A ( P I b,J| , J*1 ,36)/ 1U«B11»2H B, 

• 67*600, ,000, ,000, 65*000,999 

« >000, ,000, ,000, .000, 

- *-2*B5 r2 ■ 7~3 , -*790, ,05* 

• *00, <00, *00, ■ UD i »00 1 

• » I 0 1 53b ) 3E0 ! t - , l B50?U68E- ! , 

• -*3690ab00E'3» , l 7698932E-9 , - 


,019, ,000, 

*Q 0(1 t9?9tOSQt 

*000, ,000, 

, 00 , * 00 , 

« t 12109 3SE-2, 

. IUB 195B0E-6/ 

»C25i rO&O-j 

.000,999*000, 
•000, ,000, 

.00, .00, 
.I92Z5233E-2, 

rl 19B071-7 €-*/ 

.055, .000, 

.000,999*000, 
.000, .000, 

,00, .00, 

» r?S 599 * 9 E” 2 r — 
> 169 ) 585 1 E“ 6 / 

,1> j -ooo, 
•□or /99,ooo, 
,000, .000, 

*E'ff j - * JOi 
.281 39926t"2 , 
•29230B72E-6/ 

. 1 88 , . 000 , 
.000,999.000, 
,□00, ,000, 

, 00 , . 00 , 

• 36895955 E- 2 , 
. 3 J 819955 E- 6 / 



OQ+15- -55*— 

00115 39* 

OOJ1S 40* 

00115 91 ■ 

00115 42* 

00115 4 J» 

-00+ 1-5 9-9--' — 

00117 4S» 

00117 96* 

□0117 97* 

00117 4 8* 

>30117 99* 

00117 5 1 • 

00121 52* 

0012} 53 * 

00121 59* 

00121 56* 

— 0O+-2-J &** — 

00 1 2 1 57* 

ooi2i &a* 

00123 59* 

00123 60* 

00123 61* 

— 00+23" 

00123 63* 

00123 69* 

00123 65* 

00125 66* 

00125 67* 

00125 66-*— 

00)25 69* 

00 | 25 7Q* 

00125 71* 

00125 72* 

00|27 73* 

-60 I 2 7 79 *- 

0012/ 7&« 

00127 76* 

00127 77* 

00127 78* 

00|27 79* 

ooi3i eo* 

00131 S',* 

0013J *.2* 

00131 83* 

00(31 69-' 

00131 B5* 

00131 66*— 

00133 67* 

00133 US* 

00(33 BV* 

00 133 90* 

00133 V | • 

00(33 +2* 

00133 93* 
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- • 692592 1 ?E“2 , . 3563 l 4 24E-3 , - . 6728 1 444E-5/ 

DAHifl 38 ■ J 1 , J* | , 36 I / 0 7 , 620 , 2HsR , 1 6 . 1 05 , 1*940, ,133, 

•875, 6.863, ,000, .770, 5,563, 6,172, 
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2.68,2.73,251.300, 36.50, 

3-l*-2ttr- 2*!»4thi fr.62 , 5*12, 4»?Z. 3.59, 3*-44 T *00. 
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313* 
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it H* 
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316* 
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116* 
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00235 
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00237 
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27*300, 32 • 900 i 32«9p0i 
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■ 0 , 
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5*961 

2"»'6t7--j2*-7-3 *259*100* 

33*50. 29, '.0, 7.19, 5*58, 4,6!* 3,92, 3.73, 

rl-l77o3tvE*lt -,261317075-0, .6427 ! 033E-1 , 

-.727546S9E-2. .37467B99&-3 , - , 70793358E-S/ 

.786* 6,071* .000, .0U0, 4.867* 5,378, 

-5*-535*— 25*800 * 3Q.4O0, 3Q.4Q0, 43,800* 44 » l 00-* 0^3*000* 
*00*2.73, *DOo, 42*30, 

—46, |-5, — 26t-45t 7,79, 6*08t 5,02, 4,26* 4*03 * 

• I 1 7 9 6 608E+ 1 * -.26647782E-0, .6561 9469E-1 » 

DMAtPt 41 ,J1 *J«1 ,36)/ 92*90*, 2HN0, 18.986, 2.371 , 

-*- 7 - 46 -. — 6,724,- .000, ,000, 4,581* 5*026-, - 

5,023* 2 4 * 000 , 28,200, 28*200, 40*600, 40 . 900 * 229 • 000 * 
*OO r g* 7 3 , - — *OOo* — 45*50 * 

38*90* 2/, 90, 0*39* 6*61, 5.46, 4*65, 4*34, 4.99, 


4,71, - - 
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-.75897910E-2, . 39 I 1 1 497E-3 , -.7 39254 1 2E-57 

0«?A-lP< 42, j! rO*t ,361/ ,’5*940, 2HH0. 20,000, 2.520, ,227-* 

*709, 5 • 40 J , rOOtj .000, 4.298, 4.718, 

4-*9-i- 3-* -22*4-00* 26*300* 2«*3&G, 37*700, 38 »QGQ *z08*0QG*- 
.00,2.72, * G 5 u .* 4 8 • ’i f, •» 

•x — 9-*U6 ?»-»?*— 5*4?,' 5*05*- — 4", 68 , 5- » 2 8-, 

* * ! *r?ffc7’*»443E-0, , 6 8 2 1 3642E- 1 , 

* a*T?345334e“2» , 34 3*5 iOOE-3 , -.753568021-5/ 

OATAtfM 4J S Jf ,J B i ,36!/ 9’ * 000 , 2HTC , 2 1 . 044 , 2,677, .253, 

* ■- vi 75 1 5oiA-5* ,030, ,000, 4*060, 4*436, 

* 4*632, 2 0 9 U 0 , 24*400, 24*400, 35,000, 35*300,189*000* 

-e — ,-00'ril f7-2-v — r-QOOr'52'* - 10T ■— -- 

* 44*50i 3 1 j> 9 5 9*74 * 7*75* 6,4), 5*47* 5*03, 5*58* 

* , 118*551 iE+1 , -.28073108E-0* . 6 9 337 845E- 1 , 

* -*7B629479t i ‘2» . 4 05256 8 7E-3 , - . 76600245E-S/ 

OftTAf Pi 44, j! ,J-1 , 36 1 / I0> *070 ,2KRU,2Z. i 17, 2.83B, ,279, 

* ,643, 4*846, .coo, >000, 3,830i 4*180, 

* 4*369, — h9~* 5 0 0 ~t — 22-*8G0,- 22*800, 32*600 , - 32*-9O 0 *~l 72*0 00* 

* *00,2*72, *OOo, 55*60* 

* 47*50, 34,10, 10*47, 8*39, 6*93, 5,91, 5*39, 5*B8» 

* • 1 1 87 1 66 IE* 1 i -.285B24B4E-0, * 707 28004 E- 1 , 

* »»0i)2?2a21E-2t ,4l4|0560E-3, -,78 j06457E~5/ 

0 AT A | P ( 4b, JJ ,J«1 ,361/102. 90S, 2HRH.23. 220, 3.004, .307, 

-* *6 I 3 * 4-r69-7, ,099,- *«0O, - 3*606* - 3*44*x 

* 4* 130, 18,300) 21.300, 21 *400, 30*500* 30 » BOO , 1 57 * OQO , 

* *00,2*42, *00o* 59.3D* 

* 50*70* 36,40, 11,24, 9*05, 7,4?, 6.38, 5.78, 6.20, 

* *1 1 9 1o2-46e* l , -.29001 153E-0, , 7 1 6 1 382 2E“1 * 

* -*8l54i762E*>2, . 4205&990E-3 , -.79525062E-S/ 

J*1t 36 I/ 1 106.400, 2HpO, 24 , 350 r" 3* *-33-5-, 

,585, 4.368, ,000, .000, 3*428, 3.724, 

3r90O* 17*200, 19,900, 20*200, 28*500, 28*800*144*000* 

•00,2*72* *000, 63 ■ ID, 

50* 9 5, 36 * 70, 12,03, 9*74, 8*05, 6,87, 6*18, 6.52, 

1 1 l 927643E+1 , -.294D3836E-0 , . 7285B036E-1 , 

“ *8g 7 ' 1 1 6 78E"— 2, * 42675528E-3 , 80*94 272E-5/ 

DATAIPl h7,J» ,J>1 ,36) / 107*870, 2HAG ,25,514 , 3*351, .367, 

*559, 4.164, ,000, .Dud, 3.254, 3,614, 

3*698, 16*100, 18*700, >9*000, 26*600, 26.900,(32.000, 
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UU-Z-57 - — 

00257 367 • 

002H aae*- - 

UOZfcl 3a?» 

□ 02*1 39q* 

0U26| 391* 

Q(H4. | 193* 

00261 393* 

0QZ61 39*(* 

00263 395 » 

00263 396* 

00263 397* 

- -00263 390* 

00263 399* 

00243 - 4 00* 

00263 S0|* 

00265 40 2* 

00265 703* 

00265 905* 
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00265 H07* 
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0026-7 4+fj-* 

00267 H| |* 
00247 H l 2* 
00267 H | 3* 
00267 H1H* 
00267 H i 5* 

0027 1 6-J6* — 

00271 H 1 7* 
00271 His* 
00271 H | 9 • 
00271 H2u* 
00271 H 2 1 * 
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U0273 H23* 
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* H-20-54)686&* 4 -i -TT3253O950&-0, »-ft0986979&^l-, 

• -•92O51580E-2 • . 474606<i 1 E-3 , - , 89 6 BS973E-5/ 

-fr6-7At7M-56 fl J)-,J*l ,361/137. 340, 2HB* |37»4H| , 5,247, ,7Blt 

• *000, 2,776, ,000, ,000, 2,068, 2,204, 

- ■ --* ■—•2*562, 9,55-7* 10,845, 1 1*700, 15,500, 15.BQQ, 64*700, - 

• ,00,2,69, * QOg , 109,00, 
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D AT A | p 1 57, J| , J»l,361/130,9lO,2HLA,3B. 925, S.4B3, .832, 

-* *000-^-2^666, 14,080, .000, 1 .973, 2,103. 

• 2*258, 9,042. 10.321, ll.lOO, 14*700, 15*300, 6 q,200, 

• 97* BO, 70*20, 23*70, 20*40, 16*90. 14*23, 11*93, 10,77, 

• .12073441 E + t » -*3314545BE-0, *8260)01 l £— -1 , 

• -.93905069E-2, . 4B4 I 3 1 5 1 E-3 , 9 l 473894E-5/ 

OAT-APPI 5*,J|,J»1 ,3617140*120 1 2HCE»40.443, 5,723. ,883. 
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• -*96508844E“2f . 492226 1 0E-3 , - . 9297 1 908E-S/ 
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DAT A -fP [ — 6 4 f j-li' i}* 1 13 61/ 1 37 *-250 1 ZHG&-|50? ?39r" 7-*243 t - 1-*1B5i — 
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.00)2.66. .000)157.00) 

13-4.00. 96, 3D. 3jj , 3 1 i 29*60) 24.50, 21.20. 17.02, 14*30. 

■ 12127247 E + 1 1 -,3H91?3 2 2E-0, ,07066293^-1, 

DATA1PI 65,31 |J*1 .361 /15 b «924,2htB .51.996. 7.514. 1.241 , 

tQQO 1 - 1.9771 10.000. .000, 1.421. 1.501. 
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1 1 - 212731 4 C -* 1 ■ -v35-&-7-46-27E-~£l,- - *673 7 2 1 8 4 E — 1 1 
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vEhfr,— 32*60)— 2-7*00 .-23.40. — 18*70-, — Hr. 30. — 
. 1 2133589E+1 , - , 3530 l 3s0E-O , , B 7932 l 1 OE- 1 , 

-,99672U05E-2t .512263 745-3, - * 96536232E-57 

0 A7 6 ( P ( 67 , J 1 , J* t , 36 1 / 164 .93D,2HH0.55.618 , 8.071 , 1.351 , 
,000) 1,845) 9,200, .BOO, 1.317, 1,390) 

1.535, 5,020, 6.453, 7*128, 8.899) 9.155, 31,300, 

151*00. 108. SO, 39*76, 34*3d, 20.50, 24.50, 19.59, 15.85, 

, 12139149E+1 , -.35518014E-0, BBH59D 1 5.5- 1 » 

-« IQQ22942E-1 , . 5 1 492 1 9 7E-3 , -.970030045-5/ 
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tOQO, 1.784, 8.820, *00D, 1.260, 1 .33B, 

.00,2.84, .000.184.00, 

157.00. 113.00) 41.69, 35*90, 29.B0, 25.70, 20*50, 16.30. 

• 121460445*1 , 35747602E-0, , 990 327 1 OE" 1 , 

-.100852U3E-1 1 ,51 79590 I E — 3 , - . 97549555£-5/ 

OATMPl 69 , J 1 , J*i ,36} /160>934 ,2HTH,59.390, B.6H0, I.46B, 

- Tta UP -i 1 1 7 27, -8.490, ,000 1 t,222| - tfMfli 

1*433, 5.366, 5.931, 6*559, B,|5<;, 8,433) 27,000, 

.00,2»64, ,000, 192.50, 

165.00. 118.00, 43.69, 37*70, 31.20, 27.00, 21.44, 16. 85, 

•121488196+1) -.359279226-0, . fl 943253BE- 1 . 

I0123145E-1, .519570206-3, - . 9 7 0Q234OE-5/ 1 

— Dft J (-A-HM— 74Hd^-rO»I-T*<r) / 1 7 3*040) 2 Hy 5 ,6 1,332,- 0*444* ~Hr&38* - 
.000, 1.672, 8 , 149, >000 1 1.182, 1,243. 

1.3 66, 5.161, 5.696, 6*33|, 7*836, 8,086. 26*200. 

•00,2.63, *000,199.001 

170,00)122.00, 45.79, 39*60, 32*70, 28. 30, 22,40, 17.45, 

. t2l53540E+l , -.36129O34E-0, . S9909077E- 1 . 

■ — -*- ltll -7-203-trE— I. ,5616291 2E-3, - . 98 1 88940E-5/ ~ ■ - - 

DAT A t P 1 71, J) .J»l, 361) 174, 97D,2Hi_U. 63. 314, 9,244, 1.5B9, 
.000, (.620, 7,040, *0Q0, 1.140, 1,199, 

1,341, 4.972, 5.475, 6.112, 7.545, 7.777, 24,700, 

.00,2*63, ,000,706.00, 

1 76* OU, 126.00, 47.92, 4t»3o, 34.20. 29.60, 23.40, 18.00, 
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..02U9096E-1 1 .S23HJ868E-3. - . 9 G 4 4 1 7 33E-5/ 

DATAl Q t 72,J 1 , J"1 ,36 I / 17P • 90 ,2 HmF,65.351 , 9.561, 1 .6*2 , 

1,570, 7.539. .ODn, 1.1U0, 1.1S5, 
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0033 | 


527. 

529* 
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537. 


00333 
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592* 
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— *-*HJ22ol22E-l » .52380 10IE-3, -.98366072E-S/ 

DATA|pt 73, v |),J**l,36)/ia0,998,2HTA»67.'ll6* 9*881 , 1.735, 

> 00 (1 1 1.523* 7>2*Zl »0tm r 1*061 l — - 

1 c 255 * 9*569, 5*015* 5*637* 6*890, 7* 120 * *>4.*1Q0, 

fGB-Ttt«-62i — *000,222,00, 

190*00*136,00* 52,99, 95*20, 37*505 32.90, 25.98, 19*05. 
-»1-21&9?10E*1 * -*3662?272E"0, .9Q982919E-1,- 

-« 10i69132E-l » .52S69993E-3 , - * 9875 1 925E-5/ 

•000, 1.976, 6*703, *000, 1.025, 1.075, 

T T 2"t6-r 9*905, 9,023, 5,951, 6*690, 6*871* 20-*900, 

,00,2.61, *000, 231.00* 

*-1-97t00t192 1 00, 59,79, 97*20, 39 , [ 0 * 33.80, 26.57, 19,60, 

* I216Q669E*! » -.36780Q97E-0 ■ • 9 l 280 3 0BE- 1 * 

■ *, | Q2? 9 3 9 3E - 1 F — T52658&96E-3 , — • r9S85&79-0 

DATA(P( 75, ,361/186*200, 2HhE,71>676, 10* 535, 1 .BB3* 

.000, 1,933, 6,727, *000, *990, liQ37* 

| . 1 77 * 9*231 , 9*625, 5*253, 6*357* 6*599 , 19*800* 

.ag T 2»6t, .D00.239.00* 

209*00,197.00* 57*21, 99*30, 90*90, 35.30, 27*68, 20*10, 

»• I032S&2SE-1 , .527707995-3, - . 9 90283 I OE-5/ 

OATfttft 76, 01 ,J“l ,361/190. 200, 2«oSl73, 070, 10.871 , 1.960, 
*0D0, 1,391, 6,990, .ODD, *956, 1.001, 

1“. l“90r ' 9-* 069 , 9*993, 5*060, 6.(06, 6.357, 18*800, 

•□0,2.60, *000,297*00, 

. 1 2169927E + 1 , -.37103956E-0 , , 9 1 97 56 B9E- 1 , 

-» 1 0^85 1 09E“ 1 , .SZ891276E-3, - . 9 9 1 77905E-5/ 

DATA 1 P < 77, j) ,j«l ,36)/192.2D0 l 2HlR,76* III , 1 1 .215, 2,090, 
.000 , 1 , 38 i , 6,262, *000, .923* .967, 

1 * 106 , 3*915, 8.273, 9*873, 5*880, 6*099, 17*800, 

— ,00,2*601 ^ 0 0 3 t 25 6* DO 1 ■■■ — — — — — 

219*00*157. r - 62*27, 53*70, 99*50, 38.90, 30*01, 21*35, 

•1216998 r 1* -,37232566E-0, . 922067 1 7E“l . 

-« 10369639E-1 , *5291 6 6QBE— 3 1 - * 99 l 50B39E-S/ 

OAT A-1P1 78, dl, J** I, 361/195. 090, ZHpT, 7S. 395, 11*568, 2*122, 
*000, 1.313. 6*097, *0U0, *893, ,939, 

1 * 072 , 3 1 7 7 4 , *hrt 1 9 t 9.7-JO, -5.665,- — -5.880-t — -16* 900 1 — 


.00,2*59, *0OQ,263,00* 

225»(JQ, 161 ,00, 69,92, 56*00, 96*90, 90* 1D| 31*22, 21,95, 
. |2166262E*1 t -.37305379E-O, • 9252 t 3 1 0E- l , 

103961 12£“| , .530QU972E-3, - , 9926 1 9 15E-5/ 

DATA (PI 79, j),J**t, 361/196* 96 7, 2HAU,80*7 7.5, II* 919, 

— — . 000 , 1,27 6. 5,890. - *O0D, *863,- — ,903-,- 

1*090, 3*620, 3*939, 9*822, 5,915, 5*629, 


•00,2.59, .000,272*00, 

Z32>D0il67,o0, 67*68, 58*30, 80.30, 91.80, 32*97, 
» 12165537E+I s -.3750*9625-0, . 92727299E-1 , 

10907273E-1 . ,530l3023E-3, - , 99 1 90720E-5/ 

OAT A ( P 1 — 80-,-d )-rd-®t ,36 1 /200 .590 • 2 Hh< 5 , 83 * 102, 12.209, 
*000, 1 *291 , 5*6*l8, *0uD» *035, .872, 

t.000* 3*982, 3.779, *1*399, 5*280, 5*913, 

• 00 | 2 * 58 • • 00d ,281* DO > 


2.206, 
I 6 g QUO I 
22*50, 


2 ■ 295-. - 
15.200, 



* i u 'iw» u f 4 , # □ U3| 

00351 585* * .OQD, l.lHi .DQO, *000, .732, .7*3* 

OQibi SB**- ■* -^97, 2,957, 3.228, 3.76R, 9.993, 8.625, 12,500, 

00351 SB 7 * • .00,2,56. »OO0»3l6.OQi 

00341 508* - ♦X7tH-B0-rV9*).OO> 83.SI , 72*00. S9.70. 51. AO, 39.19, 24.30. 

f> 00351 5a?« * .J21 A2 1 82E* l . 380709 7SE-0 , .9363598 E-I . 

~tH -003-5-1 5*70* • ,&29181-6BE-3 b — ‘ t 9B 58069 9E - 5/ 

■*" 00353 59 l • 0 A T A f P 1 B5 , J ) , J* [ , 3 & ] 72 1 0 . 0OD , 2H AT , 95 . 730 , I 8 .21 9 , 2.7B7, 

00353 592-* ♦ - ,000,-1,085, .000. .OOO .709, .739, 

00353 593* * .672, 2*976, 3.099, 3*635, 8. 279, 9*959, 11.900. 

00353 599* * .00,2,55. • 00c > 32 5.00 » 

00353 595* •278*00.199,00. 37*27, 75*2D» 62.30, 53.90, 9D«69, 26.60. 

00353- -W4*- * .-l- 21S 999' t - e +l- r --~»3-9l-5»6g7£-0. r*93700707E-l-r- — — — - 

00353 597* • 10H35U35E-1 , .5230U29E-3. ~. 9B Z65 1 55E-5 7 

00355 59e« OATA (Pt B6,J> ¥ J»i ,36)7222. 0Q0,2 HrH ,98.909, 18.6 1 9 , 2,892. 

00355 59 /• • .000. 1 .057, .000. .000, .000, .715, 

00355 60u* * tB 99» 2,767, 2.937, 3.512, 9,|lBi 9.309. U.900, 

00355 60J o . , 0 D , ,00, .000, .00. 

003&5- 602^— *atrtrr ftC , 2ft r - r Z'Sr 'Hr^t r -7fi.5er 65 • 20 , S&e 30 r-+Mi7"WrWr 

0Q355 603* * , 12,5703 1E*1 . -.33295738E-0, , 9 37 8 1 7 1 0E“ l . 

00355 608* * ", )0M2866OE-1 , , 52693 1 37E-3 , 979S7573E-57 

00337 605* DAtfllPt fi7 , 0 ) , J * ] , 36 ) 7223 . 0 00 , 2HFR , 9 9 , 9 99 , 1 5 . 03 1 , 3.000, 

OU357 6U6* * *000 r 1.030, »0OO, *000, .000, .000, 

00357 607* • .825, 2.672, 2,876, 3.3B7, 3.9*0, 9.186. ID. 900. 

00357 AOU-* -« — r OO T <0Q | , OOO t .00. - - 

00337 6U9« * *00,210,00. 9 5 ,35, 82-20, 68. |0, 58.90, 93.69, 2B.30, 

00357 610* • *1 215JI98QE*! . ~.3B3 1 7627E-G , . 9 3B I 22 B7E- 1 » 

00357 6,l» • 109l93«0E-t , . 5Z592775E-3 , -. 9 7S62306E-S7 

00361 612* DAW-PC SB, 0K.0*| ,361/226. 000, 2HRA, 99, 999, 15.998, 3*105, 

00361 6 1 J * • .QuO, 1,005. .000, .000, .000. .000, 

003*1 *19*~- • *-9U3-, 2-. S/7 .-- 2.773, 3.279, 3.B26 r - 9.01-2 r- 1 0 -. 3 30.- 

00361 6(i* • .qO, .00. .000, ,00i 

00361 616* • .00. ,00, 99.70, 86*UQ, 71. 20, 6|.Sfl, 95. 20, 29.60, 


00361 617* « , 1 2152/ 10E*I » -.383B9122E-n, . 9 3 89 1 732E- 1 , 


si mid TOB^iao 


— WQi4-i— ' 
00363 
UQ3&3 
00363 
003*3 
00363 

00~3*3 — 

003*3 
00365 
00365 
00365 
00365 
- -303*6“ 
00365 
003*5 
00367 

003*7 

00367 
““000*7 — 


00367 

00367 

00367 

00371 

00371 

0037 1 
00371 
00371 
0037 1“ 
■t* 00373 

tJi B 03-73 - 

</' 00373 

00373 
00373 
00373 
00373 

— 00375-- 


619* 

* 20 * 

621* 

* 22 * 

623* 

625 * 
-62-6* 
627* 
626-* 
629 • 
■430*' 
631* 
*32* 
633* 
*3*1* 
635* 

637* 

- >30* 

639* 

*4G*- 

*41* 

6 '1 2 * 
643* 
*44* 
*45* 
*4-6* 
697* 
£»*fd*"^ 


649* 

650* 

651 • 
-642* 
653* 

w 3 , 'T * 


00375 
00375 
00275 
00375 
00375 
1303-75 
00377 
00377 
003 7 7 
00377 
003 77 
00377 


6S5* 
*56* 
657* 
650* 
659* 
■6*0* — 
66 ( ■ 
662* 
663* 
664 * 
66b* 
**** 


■ « 103 9 9 - 707C -- 1 « . 523 8 9 *896-3 , -*-9 7 1 *0 320E-5/ 


0 AT A (PI 

*000-,- .-9-80 » 

.702) 2.989. 

- tOOt — ,00, *000 

*00) .00.109 

H21 ‘ I7 - 6 - 9 -6 & * 1 t- 

-• 103*6131 E-l » 
D+T*-l-fM -“9U t j l xd-l 
.QUO) .956i 
* — .76 T, 2.901 , -■ 

aO) .oo. <oOo 

. 00) rO fl -rt-09 


, rmrr 

, 12199961E*! , 

■ — -1033 &-1-6 5 S— 1 ) 
DAT A t P l 9 1 . J ) ) J” 1 

*DOO. .933. 

•791, 2.313. 

■ a°, — *00-1 r OO fl-T 


,36 1/227. 000 »2H AC .99.999. 15. 87 1 i 3.219, 

.000, .ODD, .000) *000, 

2.672, 3*170, 3*689, 3.879, 9.917, 

V i OD * 

*30, 89*90, 79.50, *9.90, 96. BO, 30. BO, 

-w, -389-225-7 95-0-, ,9-393664 9 E “ 1 I 

•52125077E-3, -• 96528 447E-5/ 

,361/232.038 ,2HxH, 99,999, 16.300, 3.332, 
9.138, *000, *000. .000. 

2,577, 3*000, 3.573, 3,795, 9.992, 

, . 00 i 

tdtOt — 99*30, — 7-8 ,-1 0, — 6 7 , 5 0~r 98 ,9 3; — 3 1-.30 


-.3047O334E-O, . 93676 1 SOE- 1 , 

i 5 1 09O8O3E—3 , 95955593E-5/ 

,361/231 *000, 2HPA, 99. 999, 16.733. 
8.022, »OoO, .ODO, .000, - 
2.979, 2.900, 3.993, 3.619, 

-rOO-r 


3.992, 

9.098, 


*00, .00*119 

*1-2-1-9-18 67 E+T, 

-«lQ 31 o 3 a 2 £-l, 

DA-T-M-fM — 92 , j-) , 1 

e ODO I .911, 
.7- 2 -2- i 2*2 - 3 5 , - 


•50, 98*70, 61.80, 70.70, 5Q.10, 32.90, 
-.365I8625E-0, .93618155E-1, 

•516S6738E-3, 9530 1 63BE-5/ 

,361/238. 030, 2H 0,99,999,17.166, 3.552, 
3.910, .000, .000, .000, 

-2t3-98t- - 2t-&82f 3*3110“,- 


•□0, .00, .000 

--.Ott,-- rOOil2o 

•I2136379E+1, 

*t02*9*3*E— 1 r- 

OATAIPt 93,J),J=i 

— ' — rOU-e -1 .8 89 , - 

.709, 2.1*0, 

- .eo r - .BfjT -.ooo 

*00, .00,127 

* 1-21 3-272 1 E* 1 » 

I 02 1 9 3 1 8E- l , 


. OqO i <068, 

»0GO, 2.090, 

■QO, .00, .000 

—DO .00,135 
• 12129555E*! , 

■ ^ l ' OT -6 s o g s e** 1 , 


. 00 , 

,70,109*00, 86*20, 79.50, 51.80, 
-.38523977E-0, . 9 3 900300E- 1 » 

• 5 1 28 *025 E— 3 , -,94524676E-5/ 

.36 1/237. 000, 2H^P, 99. 999, 17.610, 

rOOG, *‘QOOr~ *000-, rOOO, — 

2.313. 2*790, 3.220, 3.387, 

, *00, 

*10,109.50, 90.00, 7B.5D, 53.55, 
-.3S541490E-D, . 93224943E-J , 

.50999719E-3, 9 372455 7E-5/ 

,36-17-244*000, 2HpU , 99-,999 t hB-057, 

.000, .Ooo, .000, .000, 

2,229, 2*7ib, 3.114, 3,279, 

. 00 , 

• 40, 1 16*50, 96,70, 33,60, 55.34, 
-.385S4609E-0, . 9302S923E - I , 

*5057-773-1 E— 3 , -t4286S595£-3/- 


33*60. 


3,666, 


8.264, 

39.90, 


7.996, 

36.40, 





D AT A { P l 9b, J|,J*1, 361/243. 000, 2H AH, 99, 999,18.504, 
-OOO, ,048, .000, »Uqq, ,000, .000, 

*030, ■ 000 , ,000, *000, .000, .000, 

•QOy ,00, .OUQi .00, 

, 00 i ,00, *00, .30, *00, .00, 57.16, 


3.807 , 
,□00, 
.00 , 


00377 

067* 

* -• 1U082912E-I , 

.5001 6353E-3, 9 1 6 1 0882E-5/ 


00401 

*68* 

D AT A fP t 96, J } j J* 1 

,36 1/247*000, 2HCH.99. 999, 1B.930, 

3.971, 

00401 

*<>?• 

* >000, .000, 

,ooo, .(ion, .ooo, ,oou. 


0043 1 

670* 

* *000, ,000, 

■000, *000, *000, .000, 

• 030, 

00401 

*71* 

• .03, .00, .000 

, • 00 , 


00401 

-67 3* 

* TOO i rOth 

■ 00, *00, »00« ,00, 59HJ2, 

.00, 

0Q4U1 

673* 

• • 12I20562E* 1 i 

385 1 S3 09E-O , . 92377782E-1 , 


00401 

674* 

♦ -. I0029525E-I , 

.4961 97966-3, 9O&B6890E-5/ 


009U3 

675* 

□ AT A t P t 9 7 , J ) i J* l 

, 36 1/24 7* OOO ,2HBK, 99. 99 9, 19.462 , 

4.132, 


9 S - t » 


- QB 4 U - 1 - 
004U3 

otwoa 
ooHua 

00403 

OD*toa 

_n n ti rtfi - 
UUlUa 

00HU5 
00405 
00405 
004U5- -- 
00405 
- » ■• "" • O - Q-^ SS 
00H07 
-QQ407 
00407 
0Q4Q 7 
00*107 

t}040-7 — 

00*107 
QQ41 1 
00*1 1 1 
oo4U 
001 1 1 

004 H-- 

0041 1 
0041 1 
00413 


-6' 7 ' fr * 
677* 
if - fr *- 
679* 
-6 ft 0 » 
iUl » 
-6B2’ 
603 * 
4B4* 
605* 

406*- 

407* 


■6 00 * ■ 
609* 
6 * 0 * 
491* 

* 90 * 

693* 


_ 4^x| t 

69 s * 

*96* 

677* 

696* 

699* 


■700*- 

701* 

? 02 * 

703 . 


«OUOi .QUO, 

* — r£lQ| — tOO, *00Q 

• DO * . 00 , 
-»t2Hr&70E-*“l T 

• 9 9 b 0 6 4 96 E“2 * 
A - tlH - 

• 000 i • 000 i 

— rOOOl — » OQOt 

.00, .001 *000 

rOQr ,OGt ' 

. 121 16747E+1 • 


>■ 008 -1 


* OOO - t " 




DATA { P I 99,Jl,J"l 

.-QQOnr *000* 

. 000 • .000, 

— *-0Qy — .00 » -—*000 

* 00 , . 00 i 


r-t- 2-l - l - 2 ^'9-4- E *t-r- 

-.97807777E-Z, 
i>AT"A iP-H OQ* j ), j» 1 
> 000 , . 000 , 

’■-*■000-1 iOOQi 

• 00 , . 00 , .000 


* 000 * * 080 , 

.000, .000, 

, « 00 , 

•OQi *00, *00* 

’-.3B51-69I5E-0, 
.49175555E-3 , - 

0*2 

,000, * 800 , 
.000, *000, 
i * 00 , 

* 00 , ■ • OQ , • oO » 

-,384Bt9e7t~Q, 

* 40 - M- Q4-6 J i&-3*- - 

. 361/254. OOO.ZHeS. 

.OOP, *000, 
•000, .000, 
v . OD , 

.00, *00, *00, 

—*304-1-64435-- 0* 

.H7936640E-3 , - 

,36>/253*000,2HFM, 

.000. >000, 

- .000, *000, 

00 , 


* 006 , * 300 * — 

*000, *000, .OOOl 

*00, 60*93, *00, 

• 92H3495E-1 » 
.B9659S07E-5/ 

• 000 , > 000 , 

*000, -»000* — *QQOr 


, 00 , 62 * 8 0 *, 

,91711 203E-1 , 


■* 00 , 


99 i 9 99 ,20*410, 4.37H, 
• DOO, .000-* - - — 

. 000 . •□□□, . 000 , 


• 00 i 


> OQ i 


— rOO-i r 

*1211 1 34QE+ 1 * 
-*969201 7 0E~2. 
END 


— rOOT - 


- o-OOl 


• trG , 


.00. 64.87, 

. 711 8 0729 E " 1 , 

.86B73597E-0/ 

99.999, 20*900, 4*498, 

.000, *000, 

.0001 *000, .000, 

rOO* 


-.3Q3636B3E-0, 
.H726906BE-3f * 


rQS-j — 66* 70 | 

,90 6 93355E” 1 , 
•8537Q94BE-5/ 


end Of couPIUation: 

- DBtOC-K & YHBOfc - l - e —' — 

oauocn code keuocma&le 


NO diagnostics. 


34-30 Q — 0+&75535- 

19 NOV 71 10;2D130 1 0J6H1227 

0 01620707 


1-4 ■ 
1 H 

m 


703 (DEt-eTElH 
1 IDEUETED) 
600 


UJL-NOV-71- 




H+S4T-W-S-! 

UNlVAC uoa FORTRAN V LEVEL 2206 DU24fl -IEXEC0 LEVEL E12Q10009A1 
THIS C(jhPlLA--fTGN— 4-A-S DONE— ON — l- 9 NOV 71 AT 10JZU14 


SUBROUTINE K15J ENTRY POINT O0O271 


STORAGE -USE 04— EODElH -BUD.34-2..; -DAT* (Ol 0004601 BLANK C0HM0NI21 000000- 


EXTERNAL REFERENCES (BLOCK, NAME) 


COO- GRDPI 

oooR— plot 3 — - 

DOBS UHPBUF 
BQO*- - N*DUS 

00(j7 N 1 0 1 3 

OCUl N ERR 3 s 


STB RACE— ASS IonhENT <-BlOCK, TYPE* RELATIVE- LOCATION, NAHEl 


~"~*OQtTf 

-I-2-IG- 

eeo-t 000054 - 

133 G 

ooo-i- 

— 

-000-132 15AG 

DO 0-1 OB 0-1 5 o — l-65-fr — 

0SO- 1 DGq-J-A-J— 

-142G- 

good 000337 

20F 

OOUl 00020 7 

200G 

0000 


000375 2 | F 

Oogl 00023s ztis 

0000 00 Q 4 I 0 

23 f 

QQ{jQ -OP 033-3—3 OF 

— QOOJ 000064 

4L 

GtlO I 


□00 1 0 1 SL 

0000 R OOOQDO B 

aooo r oooi^s 

Ff 

DOqO R C 10P31S 

FMT 

OoUO I 000322 

1 

DOflO 

1 

000332 1NU 

0000 000H30 INdPS 

0000 I 000312 

IR 

0Ot]U-- i -0P033 T- 

ts- - 

0000 “I -00032 I 

J 

DOOD 

I 

00032? NN 

0000- r 000323 NT- 

ODOO H 000326 

R 

DOuU R 0Q032R 

STEP 

OOUO R 00D 325 

VK 

0000 

R 

00033D XX 





tn 

■*4 

□ □JUl 
ooiui 

— - !■*- 

2* 


■■■■ 3*' 

OUIUI 

H» 

aoioi 

5 3 

00 1 U3 

6* 

90 1 03 

7 , 

oaioH 

8* 

001-04 

•*+ 

001 12 

10* 

00112 

m- 

001 U 

1 2* 

QDl 12 

13-* 

001 12 

14* 

00114 — 

~ Jb«r — 

OU 115 

16* 

PD l l 7 

17- 

00125 

i a* 

00 1 26 

IV* 

00127 

20* 

00127 

— 2r» — 

00127 

22. 

00127 

23* 

00 1 30 

24* 


SUffRouTlNt HI SI (X»N,NG*XH!N.XRAX,F, INXiTJTLE) 


-c CALCUL ' A T E S — AR C- P LOTS 0 >r THE- Pftl-NTER— A~ ttTSTOSRA ^I F R O M— A-S-&T 

C OF DA~A POJNTS FROr GROUPED OR UNGROUPED DATA POINTS. 

0 I HE NS I ON X{R](F|NG I * B ( 1 D I I »FF I 1 D I > ■ I R I 3 > , FHT I 4 ] 

t i title im 

DATA FHT t U/6HI I 3X I l/lFHT I 2> /BHPEa*3i/tFHTl3)/ , 6H2X /,FHT(r)/6H 

i A-rt-fri OTTY ItJl *-lH !-/ - _ 

DATA ( 1 R I I ) |I«1 ,3|/0, I . 1/ 

WHWWHrwHrHM «■«* mt m »**•*■*■* m m m b* «+*«■ «.**•*« •* **«•■*«■ b wt 

C TESTS AHETHeK THE DATA POINTS ARE GROUPED OR UNGROUpED. IF THE LATTER* 
t CALLS’^GROUr SUBROUTINE TO GROUP THE DATA POINTS. 


— -WT— ’ ■ — IT — ■ — 

IF IlNX.EQ.O) NT - & 

*R-HE-tHT j30» TITLE 
3U FOrHaTI 1HI * / / / * 5 0 A * 12A6|///I 

call gropita,n»ng,xmin*xmaxiM 

STEP ■ UMAX ~ XHINI /NG 
FINDS rtAXlHUH FREQUENCY A«0 COMPUTES THE K FACTOR F0« EXPANDING OF PLOT. 



VK«1 



03 


00 1 32 
UQI-35 
00 1 37 
-00-1 44- 
OU 1 42 

-oa 1 43~ 


00143 

00(43 

00143 
00 1 44 

00144 
-004-44- 


-26*— 

2** 

27*- 

2a* 

29* 
30* 
-3-fc-» — 
32* 
--35*— 
34* 
35* 
36* 


— ft-* mi 

DO 3 I *2 | NO 
4 'fi*/r4/rX-HSrf Mil 


4 IFt<R/N*VM.GT,l,} GO to 5 

VK«VK-.|. 

GO TO 4 

-6 -V 4 -«r V ie- H 


e — vtrK-i 1 ou-t— tftfiuetieNTS and- heading for plot* — -— 

c— — — — ------ — 

20-FOfi H ArT-l-20*7 5H6-T E P ■ p 1 5 » 8 / 2 0 X , 2 9 H C E H T E R P 0 1 N T- 0 F -- I N 1 1 4 * tr~6 R OUP ■ f 1 5* 6 T 

1 I OX i 2UHN0 • OF OBSERVATIONS.^ /20X , 27HCENTERP0 J NT OF FINAL GROUP 

— 2*f 1 5 » B ■ 1 0X -- . - 1 4 UH0 I 0? — aft041P6- (6 t>20Xt?H K - r*eT Oft* r* ■ 0 // - 1 


0D145 30* WSlTElHT^OJSTeP.AHIN.N.XHAX.NG.vK 

QOlbS- 39* DO— 1 — t-*4- r ttn - - 

00 1 AO 40* I FF [ 1 1 ■ I “I 

00 l * z 4(-* 24— f4)R MTT-t-^A-rl OHPERcfcNT AGE , SX T 1 Hu-r | Q 1 7 X ■ I 31/23 Xr 1 04 l-rt*r9H********* 4*V 

00 1*2 4 2* 1 H* 1 

-004-41 46* #ft ME 1 HT r 2 - H— HtI - 1 0 rtW rHll — — — 

44* C--— ———————————— — — 

-45* C— cO«PUTE5--VAtOEFOR EACH GROUP AND PLOTS THREE LI N&S- FOR EACH- GROUP - *1-TH 

46* C THE VALUE PRINTED OUT ON THE SECOND LINE. 

40* 00 2 I f J i NG 

49* Nft*+0 O * VK -*- t - F Ml / N 1 


001*3 

001*3- 

001*3 

004*3 

00171 

-00-1-74- 


0017b 
00176- 
00177 
0 0202 
00203 
00203- 
00203 
00203 
0020 * 
QQ2UA 
Q02U7 
— 003(0- - 
00217 
00220 


bO* 

51* 

52* 

63* 

54* 

35* 


XX*XH t N* ( | -M *STEP 
ISwQ ■■■ 

DO 2 j a i j J 
!Np*(fH il > 

2 CAM. PL0T3MND,IS»XX.a»FF,NNiFHT,INXl 


5** 

57* 

SB* 

59* 

* 0 * 

*2* 
b 3* 


C WRITES OUT TRAILER INFORMATION ON GRAPH. 

23 FORMAT 1233 > 10MH*,9H*********1 , IH + /8X, 1 OKPERCENT A<»E, SX >1 HO . IQ<7X, I 
13(1 - 

«KlTEl^T,23) U.l*10,ia0*l0l 

( F | H 7 tC G » 17) C-A-ffc — P4PBUF-~ - — 

RETURN 

END- 


END Op COMPItAT-I-Ottl 

Hist SYHBQLK 

n- 1 -s L trout Rf- mc ATA &rfc- 


ND DIAGNOSTICS. 


19 NOV 


71 

-H- 


10;2D:3* 

-lt((20J-36 


01*44273 
-0 J-54-7-1-17 ■ 


14 
— 2n 

14 


*3 [QElETEOJ 

1 tDEtrrcDt- 

39 


U-HP & , 1G1 V E 

UNIVAC iJQS FORTRAN V LEVEL 2206 a02**A -CEXECB LEVEL E 1 20 1 0009 A ) 
THIS CoHPltA-TION— #A-S-OGNE— ON— J-9-N0-V— 7 1- AT 10521! 1* 


- -19 NOAt 7-L_ 


— ltW-2+H-6 T 4«7- 


SUBR00T1NE l&IVE ENTRY POINT 000353 


ST0RA6E-- US60 f-coOE-H-) — GQDS-2A-1--DAT A t U ) QQ0J*5j BLANK CDHHONtZT OOOOQU 


EXTERNAL REFERENCES IflLQCK, NAME) 


0003 CROP I 

quqr ptcria 

OOoS OHPflUF 

-000* FrAinrar — 

0007 N102S 

-ou-j-o- — -tri-tn-s 

out l N£KR3s 


S-T-OH AS E— AS S-ifiNM&NT- — t-BtOSS-r T-TPE-r- RELATIVE LOCATION, NAME* - 

— — — 80j}0 — — J-OF -OOQQ- — -GG0&-33— ITH-F OOOl GGG036— T2-*G 000-1 QBG 07 fe ~H -6G OQOT OODiOS- TS2€— 

OBOl DDOUf 1*7G OOUQ OQOOH6 20F OOOl 000212 20 I G OOoD OOOlOo 2IF ODDI D0Q2*»3 2J2G 

QOol - 0003-1-3— 225G- OOUO — 0DDH1 23 f GOD! 00D32Q-2310 - -OOOl 00012** fiL OOOl 0001-3-l-*L 

OOoU R OOaoOQ 0 0000 R aOQOOl FF 0000 R 00D01R FNT 0000 J 000022 I OOOO I 000002 IF 

-0Q00--I ODOoaa-TND OOOO 00 01-3-i - J-NJP* 0000 - 1- 000 1-30— lft- 0000 -J 00002* 1 S 0000 I 00002** 1 1 

^ OOOl! i 000025 1 2 0003 I 000U3Z J 0000 t 000027 NG 1 0000 1 000023 NN 0000 I OOpnZO NT 

■or— -0 P 0 H— 000-1- 3 o— R 0OOfJ-ft-O BQU-M— SteP— 60 0 0-R-00003-1— X-X — — ~ 


OQIOI 1* S^OUTlNE 1G1 VEU,N|NG,XNlN,AHAX,FitCUH, INX> 

■ 0 0- HH — 2r*~ 


BB10I 3* C CALCULATES AND PLOTS ON THE PRINTER A CUMULATIVE FREQUENCY POLYGON 

- 00101 -N* C — FROM— A'-SE-T— OF D-A-TA POTMTS, — - 

00101 5« C-— ———————— 

OQIUA 6*- 0 THEWS I Ow X [ IT ,F ( 1 T , 9 u 1 »FF III 1 1 F «1 01 i FHT HI 

0010** 7 • EBUlVALENCE (R,lR) 

0010** 9 * c FILLS The B ARRAY #jTh *x* to BE PRINTED for the plot, 

BO l 05 11* DATA FMT t 1 I/6HI 13X , [/,FMT I 2 >/6HPEB.3,/,FHT 131/6H2X /»FHT*«t»/*H 

00105 12* - 1 / t B t i ! / 1 H X / 

00113 13* NT* 1 7 

“00-H-4 - — “H* PF-t l -N- X - rEO -,- 0 -)— NT** — -— — — 

OOli* IS* AKITE (NT • 1 0 1 1 

00120 16* rOI- FORMAT < 1-Hll .... . . 


00121 17* #R l T£ 1 6 > 1 0 1 

00123 la* 10 FORMAT- lA//t50X,*SONCENTRAT10N CUMULATIVE DISTRIBUTION**///! 

00123 19* c — — - — — 

- 00123 — 20* C— C 0 H P UT 6 — p R &3 U irtf C-Y—Y-* t tlES — AT—EACM PERCEHT1 LE-POR-Pfii-NTQur-A T-TOP 

00 I 23 21 * c of plot 

00123 22* C- — 

0012H 23* STEP-tXMAX-XHlNJ/NG 






00130 

00)30 

QQ13U 

00)30— 

0Q130 


00)32 

00132 
00)32 
-00(33- 

00133 
-00-1-33 — 

00 1 3H 
-00)43— 
00143 
00)43- 
00)43 


00144 

00)56- 


3 I F ( t ) >1 *N/ 1 0 


TESTS #HETHEK DATA FOl«TS A«E GROUPED OR UNGPOUPEO . | F THE LATTER i 
“C-AbtS-fiROUP— SUBROUTINE TO-'BROUP THE data POINTS* 


00157 
OOF* 1 
001*2 

00-f43 

45* 

--4 6* ~ 
47* 

— .am 

-- 

001**1 

x a w 

49* 


00165 

— 50* 


001** 

S J • 


00171 

52* 

- — 

-p. 0017) 

53* 

C- 

Tjr — OI)-}-7-l 

— 54* 

— c- 

*=> 00171 

55* 

c 

00)7-1 

-56* 

c- 

00173 

57* 


00 174 

&a* 


00175 

59. 


■ — * —00 

— 40* — 

— 

00200 

*1* 



00204 

00205 
U0207 

- rtn rt.i_L 

feja 

64* 

65* 

0021 H 

6 7« 

002)5 

6B* 

002)6 

69* 

00217 

7U* 

0021 7 

7 1 * 

■-7-2*-- 

-—-•002-17 

00217 

73* 

00222 

74* 

00222 

75* 

00223 

76* 

00235 

77* 

- - 00237 

78*- 

0UJ40 

79c 


-Afll-T&S — OUT — ARGUMENTS— AND HEADING FOR THE PLOT » — — - 

-20-F-ORHA-H-20-AT5HS-TEp»pIS.S/20K,29HCENTeRpolNT OF 1 M *) *L - GRoUP-f 15»S T 
1 1 OX 1 20HNQ * OF OBSErVATIONS-15 /20X , 27HCENTERP0 1 NT OF FINAL GROUP 

WRITE l NT * 20 1 STEP , XM I N , N , XHAX , NG 

-2-*-PORN*T-(-B-X-rl DHPERCENT AGE r&X s I HU 1 1 0 t 7X r I 3 ) //fl X . 9HFRE8UENC Y ,6X|)H0rlB 

I ( 7 A i 1 3 ) 1 

COMPUTES ACCUMULATED FREQUENCIES BOTTOM UP* OR TOP DOWN 

WRITE <KTi2l> Hi j»IO*lOQ, 10) , [ I F ( I I i 1-1,10) 

NN-) - 

IFUCUH.EU.-H GO TO 5 

-- - - 

12-NQ 


5 U-NG-1 

2 z* | •- - . . 

6 DO 1 1*11 i 12, ICUH 

- l l CUH I 


connecting points. 


GROUP ANO-r toT S- EAC-H ' ) 


IR-4 

I 5-1 

ngi-ngh 

1 P H CUH . C B . | ) - N G-t-»NS -- -- 

00 2 

— i ND-0 

XX * XHJN * ! • 5 t l P 

IF rlCUN.EU*-! ) XA * XHIN ♦ II-H»ST£P 
IFH.Etl.NOl) 1R-0 

^ o-2 - j -o-rt-R — — - — - 

FF( I 1«I tS.O-Jl *FU >+J»F( Hi J 1/5.0 

- FfiH-PFt 1 )/N*iO£)» 

CAL> PLOTJl INO, tS t XX ,B,FF»NHiFHT, 1NX) 

2 | NOT' 


at— bottom— of — gh-a^-h — 


23 FORHfrT-LBXiTHFHEOutNCTiiX* 1HD. 10(7X.I3iy/*X*t0HPERCENTAGEt5X,iH0 , V 
I (7Xi 13) I 

Art I Tf- I-WT 1 23 T ' Up III i !*l * IOI * II, I ■ 1 0 , 1U0 , 1 0 V 
IFtNT.EB.l7) CALL 0HPBUF 

return — • - — - 

END 


END OF COMPILATION) 


1.0 DIAGNOSTICS. 


ORIGINAL PAGE B 
OF POOR QUAUTH 


SLi.OC*TABLE 


14-i-ve 

IQJVE COUE 






— NOV-7 | . U|4| H9 -.-ABA 

UN l VAC l IQG FORTRAN V LEVEL 22U6 002RA -(EXECa LEVEL E|2QlOQO?A} 

— THt-S-CoMP-lLA-T-lON— ^AS-DONE— UN — |-9— NO-V-7 1- AT 10121; 19 — ----- - - 


SUBROUTINE GRqPI ENTRY POINT 0{JU|H 


,sroRAfi£_usE{>4-eoBE-U-j— ooo13-J-t-oata-iqi qooo27; Blank comhonI 2) -oooooo - 


external references < block i name) 

0Dq3 NERR31 


iT^KAftt. -A^Si-fiNHE-NT— i-BLOCKt—TYFE^-RELaTIVE LOCATION* NAHEl - - - 

0UD4 0000-2-Z— BS9-J 00G0-26--H-3G-- — DO 01 DQQTOo-*fc OOQfl — j — DBOGG-I — I ODOO 0000-10— TN^F* 

flOgO l 0QQQU3 K 0[)0Q R DODOGO STEP 0000 R 0QD0O2 Z ODOO R OGODOR ZK 


O' 

M 


— — CBTO-1— 
Q0103 

00 i OH 

001 OS 
OOt 10 

^ QG 1 1 2 
- 00112- 
001)2 
00 1 U 
GO ) IS 
00 HA 
00117 
-OO-H-7— 
0OU1 
00122 
0012*1 
00l2t 
00 1 2A 
-00-1-2A- 
DO 1 2 A 
00130 
0013] 
00133 
GO I 3H 


1* SUfr»OU'T-mt sfl OP-{-Ut NrNe-*'XHIK-, Kh**^p-1 

2* DIMENSION SlUiFU) GROUP 

3a- - STEPMXHaX-XHIN>/NG 

*t* 00 b I»1,UG 

■-■■5*- - ~*-FtjUG»0 — - — - -GROUP 

4* DO 9 I*1,N GROUP 

B* C IN sTHICH GROUP Dots THE POiNt FALL 3R0UP 

IBP Z**IXm-XhINI/STEP 

■ 1 !• K»Z- - 

•DIAGNOSTIC' THE TEST FOR EQUALITY BETWEEN NON-INTEGERS HAY NOT BE MEANINGFUL. 

t-2a -TPI-X-t M r g. -0-rXWi-tt-l- K-a-j -• — - 

13» ZK«K 

H - • IFlTflfTtf ril K-K-H 

ISa IF(K*lE»0) GO TO 9 

l **■ 1 PI-KtSTtN G ) GO- To V 

17" C- — 

£ C-AL-C-UtA J f-rON— OP-FREQUENCIES- FO r-E aCH-P CM NT-FoR-E-AO«-fiftOOP GROUP’ 

19* C — -■ ■*- 

20 * — - F-iRlap’l K-)’»-l wO 

21* 9 CONTINUE 

22* RETURN ■ - - GROUP 

23a ENO GROUP 



ENO 

Of compilation; 

1 DIAGNOSTICS* 









GftoPt 

-- SYHBOfIC- 



- 1 9- NO v— 71- 

10;20i*f! 

•0 

01 AS230I 

14 

23 

t deleted 1 

SKoPi 

CODE RELOCATABLE 



i> *** 0 V 71 

io; 2 o:ri 

1 

01653217 

1R 

1 

IDELETEO) 


. 






o 

01 &53003 

»*) 

10 






B-F-OB-i » ?L 0T-3-irtHrti-F3 — ~ — — • - !-?-■ NOV- 71' - - - 1072432 1 t-546 — 

DN|VAC I 106 FORTRAN V LEVEL 2206 002*)* -1EXECS LEVEL E12CUQ009A) 

THIS £oHF|LATjUN--|lrA5--DONE-QN- 19 NOV- 7t *T 10X21 SZ1 


SUBROUTINE PLOT 3 ENTRY POINT Q0U2HD 


-S-TOfiASE-U5EO"4-COBSH-T-00026-34-DA-TA(Oj 00UZ*3j Blank COMMON<2> OOQOOO - 


external references (block, nahei 


0QD3 NWpUS 

DO OR NJOIS- 

OOoS NI02S 
OUL*- NE-RRSr 


storage assignment (Block, type, relative location, name! 


OUUl 00021 1 |0L ODUD 000163 JOIF 

- •— ouol QOQaA-tr-I-frBft-— OOCn- QOD076 162 « 

OOQl 000202 2l*iQ Ouu I OOQ2ZO 224G 

— 8BS-I 00 01 * 2 - 9 b OpOO— ft— OO p O OtJ— A 

oaoo 1 0001*1 1 0000 000172 INjP* 

OOqB R- 0001&*-S-T-AR - — QOUQ R- 000157 ZERO 


0001 

0001 

0001 


OOOOIS 1220 
000101 1-460- 
000073 


0000 I 000162 J 


pool 

-0001 

0001 


000036 1 31G 
0001-12 1 74G 
000107 4L 


OOUO l 00016Q NT 


0001 0000*1 1*I7G 

-000 1- 00q0&2 2L 

0001 00016S 7L 

-OOO t -ftS— PHT — 
oooo R aooiss plus 


■or- 

Ol 


00)01 

00103 

00104 
00107 

□ DUO 
oart-6- 


00117 

00121 

00)24 

00)26 

00127 

- - 00 13 )- 

00133 

ooi i* 

00140 

00141 
00143 

— 00444- 
00146 
00151 
00 1 52 
00155 
00160 
00 U I 
00164 
0016b 
00170 


1 • 

£•■ 

3* 

4* 

5* 


7 • 

■ 8* 

9* 

to* 

11* 

- 12 — 


13* 

14* 

IS* 

16* 

17* 

-ta- 


li* 

20- 

21* 

22 * 

23* 

■2H*— 


25* 
26* 
27 • 


SUBROUTINE PLOT 3 ( IN0,15,X*B,F,N, F RH,tNXl 

• DIMENSION A HOI 1 ,011} ,F( l» ,FRH(H1 ,FHT16l 

OAT* FHT ( & ) /6H»10lAl/tFMT(6l /JhI/ 

-fcl- P0RMAT423 X ,101X1} 

DATA BLANK/IH /.AXI5/1H1/.PL0S/1H*/,STAR/1H*/,ZER0/IH07 
NT » | 7 


-FLOTI 

PL0T1 


NT 


11 


IF ( lNX»EQ*0] 

Do- IT" j-lrv 
FHT I 1 J-Ffthl I 1 
A( l-Ofl-AXlS 

IF ( IND.Etl, 1 J A ( 1 0 1 (-PLUS 

-l-f-t'lS- r E O i- l 'l- q -Q-Fo— 2 


DO I 1*2,100 
Ai H-*BL-AnK 
A t 1 ) -AXIS 
WiKiNEtfT n 
GO To 6 


plot i 
ploy 1 

-PLOT'!’ 


At l ) -PLUS 


DU 3 |-1 ,V( , 

A t 1 1—A X- 1-5 

DO 3 J*l,9 
* 1*1—0 1-STAR 
GO TO 6 
-D0 - 5 - 1 - I , ' V -ti 


10 


10- 


A ( 1 1 -PLUS 
DOS J-t ,9 
5 A t I ♦ j ) -B lank 


plot i 
plot 1 
plot 1 

PLOT! 

PLOT! 

-ptort- 

PLOTI 

plot i 

PLOT 1 

plot 1 
plot 1 
-PfcOTt 
PLOT! 
PLOT! 

plot 1 


-.i. ■*'- 






-tfe-l+J— 28* — 8-00-7 - - - 

oat78 29 * j-rm + i.ti 

EJD177 3U* IpH JirtT-w 1 1 » OR* I «J »Gt» (01 1 I 00 to 7 

0020 J 3J» JK ( A | J) .Ut.BLANK ,0R, A l JJ *GT(QLAN|C) GOTO*? 

00303 32* 

00208 33* GO TO 7 

flOSOS 38* » "A<oK;ERtl — 

00208 35* 7 CONTINUE 

aci2io 36* I e-c-j ffO »N ?■*-*- 1 GO TO JO 

gdzu 37 * *aiTn»T.loi> ia u i , 1 * 1 , 1 □ i ) 

00220 38* RETTURfr — — - '■ 

0DZ2J 39* JU ArfVITElNTifHTfr X ■ { A 1 1 ) • I *1 1 1 0 1 > 

— — -002-30 HO* — ftfcTORH — 

00231 Hi* END 


AfcflT-h- 

PLOT! 

PLOT! 

PLC?! 

p«.on 
PLOT 1 

— — Ptt^M 

PLOT 1 
PLOT! 

PLOT! 

PfcOT-1 — 

PL0T1 


ENU Of COMPILATION: no DIAGNOSTICS* 

'PLOT 3 5VHB0M-C ‘ - 

PL0T3 CODE RELOCATABLE 


19 NOV 71 
19 NOV 71 


10:20:83 d 01*53235 u hi 

10J20:h3 I 0 1 6587*3 1? 1 

g — &j-fr5H333 18— -*0- 


tDELETEO) 

IOELETEDI 



1-?_N0V ,?4- 


a FOB,* Pc oRRE.prnuRE __ 

UN] VAC HUB FORTRAN V LEVEL 220* 0D24A -(EXECS LEVEL E12QI0007AJ 
TH]-S^CuHfi.JtA-FH)N-VHLS-t>ONE-ON— 1-4-NOV— 7|— AT 101 * 112 * 


function pcorre entrt point 0002*7 


— S-T0H-ME~U5£P4-eG0E-(-t3— 0002^*4— DA-T-A-fO T - 0O00S24 -BLANK COMMON 1 2 } -000000- 



external references (block, nakei 


0003 GAMMA 

- 00q4 DS3HT 

00 0 S DATAN 
— OOQ*- HERR-St 



STORASe ASSIGNMENT |BlOCK, TVPE, RELATIVE LOCATION, HAKE) 


UOU305 1-Jl- 
Odfn +*~4lL- 
□00007 FI 


aoou o OODOOl R2 



0001 000052 134G 

-000-1 000202 5 IL 

ODOQ D ODOOll FN(JM 

OpOQ D 000005 SUM 


□ 001 
0001 
0000 


000007 ISL 
-0002-4-7— 57tr— 
000020 FREE 


ODO) 000214 1 & Ifi 

-0001 0002-54- *0L 

0003 H 000000 tJAMHA 


OOOO D 000003 TERM 


□Opl 00004 A 31 l 
- 0000-0— 00001-3 DEMOtl- 
0000 1 000021 1 


A 0S+04 1* FONC-M-ON-FeoRR E— ( -H-rNP-FS-l 

0 *' -'03 2* DOUBLE PRECISION R2, TERM, SUM > FI i FNUH, DENOM 

— 00103 3* c- — — 

001(13 4* C PURPOSE 

00 1 03 &•-- C £«Si!rUATE-PR03 ABILITY FOR no correlation between T#0 VARIABLES 

QDI03 *• c 

— — Q04-03 e — — — — — — — 

QQID3 8* c RESULT » PCORRE IRi NPT5 1 

oai03 - ■ »• -e 

aoioj 10* c description of parameters 

00103 II* £ ft “-LrJ NEAR CORRELATION coefficient 

00 1 03 1 2* C NPTS - NUMBER OF DATA POINTS 

— .. — 00 1 03 — j-3*- £ 

00 l 03 14* C 

00103 — is* — c evaluate number of degrees of freedom 

00103 !*• c 

00104 17*- NFrEE— “-NP-TS - 2 

OOJOS 1 8* IF (NFHEeI 13 i 1 3 , 1 5 

.. 031-1-0 HM — ■ - T 3-PfrOSR C - »- tht - - - - — — 

QQIll 2U* GO TO *0 

00112 21* 1 5- R2-«-«-*-*2 - 

, 001 13 22* IF (1.-R21 I3il3,17 

DOJI*- - 23* -1-7 -NEYEN-*-2*i NFHEE/2 l 

, 00117 24* IF (NFREt - NEvEpl 21,21t4i 

e 00-1-13 25* C“ - - - ■ - — 

00117 2i* C NUMBER OF DEGREES OF FREEDOM |S EVEN 

00117 27* C 

00 1 22 28* 21 1 M a X * (NFREE-2J/2 



99H 


OO+ii- 
0012 4 
00125 
0012A 
qq i 3 1 
0DU2 
-Q5H+— 


- 27 *- 




0013* 

0 Q i i 7 
00140 
0014+ 
QQ142 
-00+44- 
OQJHb 
0014*- 
0014* 
0014* 
0014* 
-05+47— 


OQJSU 

00-141 

00152 

----- OO+fc S-- 
00 1 b* 

Q OI^.. , 

001*0 
001*3 
UO 1 *4 
001*5 
001** 
-00+7-0 
00171 
OQ 1 72 


3u* 
31* 
32* 
33* 
34* 
- 35 * - 
36* 

3 7* 
3a* 
37* 
40* 
-4 + «- 
42* 
-4-3* 
44* 
45*- 
4 1* 

-47-* 

4H* 

44*- 

50* 

- 5-1-*- 
52* 


TERM m ABS (R) 

SUM -*- TERM 
IF t | MAX ] 60 r 26 >3 1 
ii-PtORRE =— +» — TERM 
00 TO *0 

- ST 00 —3 ^— 1— 1 -r t -M- A - X — 

FI » I 

- FRrr*-* — iH-frx — J— * ■ i 

DENOM ■ 2*1 + 1 

TER M m "T&RM-* It2 I FNUM7FI 

36 SUM > SUM + TERH/DENOM 

-- f5*H«*+l Fr EE-H-i 
PCoRRE ■ 1» - PCoRRE«SUM 
■50 -To™60 

—-4HM86R-OF 0&SSEES OF FREEDOM JS ODD 




4-1- - T - nn - J t-*— f 

TERM * AbS|R1 • 0SQRTU.-R2J 
-SUM- * 0*-T-Ak | R-2 7 T t R M ) 

IF UMAX) 57 » 45 i5 t 
-SUM * -SUH— * TERN 
GO TO 57 

-5 U M-« - -5 U M- *-T E WH 
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00 0 5 DEXP 

00Q4 NERR^s- — _____ 


storage assignment (Slock, type, relative location, name) 
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quantitative analysis of RBC set 

i. INPUT-OUTPUT TABLE 
2? concentration frequency DISTRIBUTION 
3. CONCENTRATION CUNULATIYE distribution 
* t* concentration correlation analysis 
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K IS THE ELEMENT ANALYZED 
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NUMBER 

1 

Z 

3 

4 
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A 

7 

B 

9 

|0 

It 

12 

13, 

J1 

15 

1 6 

17 

18 

19 

20 
21 
22 
23 
2* 
25 
2* 

27 

28 

29 

30 

31 

32 

33 

34 

35 


INTENSITY 

COUNT 

*48, 

177. 

1*1* 

5*3. 

383. 

7*2. 

357. 

217, 

• 11 . 
413 , 
980, 
23|7, 
18*3 • 
1810 . 

1 721 * 
2*9 , 
213. 
215. 
1758. 
1«1S, 
2113. 
1*17, 
13*1. 
1126. 
18*. 
1 * 1 . 
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*00 • 
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*99, 
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,99**9 
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,99**6 

.99*70 

.99*70 

•99**6 

.99**7 

.99*72 

.99*74 
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.99*71 
.99**6 
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• 99**9 
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.99*79 
.99*«0 
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.99671 
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• 99**0 
.99*75 
.99*75 
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.99*72 
.99*73 
,99*94 
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.99323 
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.99335 
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.99339 
.99333 
.99331 
.99322 
.99332 
.99339 
•99333 
,9933* 
.99335 
,99327 
,9932* 
.99313 
.99332 
.99334 
,99330 
.99333 
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•50538— U* 
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,63822-07 
.*4340-07 
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.47*74-06 
•S*2u7— 0* 
,432*1-0* 
,3*272-0* 
.29576-0* 
,127*4-0* 
,12109-0* 
.22930-0* 
,15759-0* 
,25646-0* 
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,18199-0* 
, 15196-0* 
,71439-07 
.142*1-0* 
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• 67 
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.13 
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.07 
,34 
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.13 




CONCERN FREQUENCY D1STRJ8UTJON 




STEP* .117713*5 

CENTERPQJNT OF INtTlAU GROUP* .Q72510U NO. OF OBSERVATIONS* 35 

CERTERPOINT OP F|HAU GROUP* IUQA71572 NO* OF GROUPS* 7 

K FACTOR* 3, 


percentage 

7,251*0? 

2,203*01 

3,480-01 

S,15«"QI 
8*435*01 
a, 1 13*01 
9,590*01 

PERCENTAGE 


0 10 

t- IP «»•*••• f + 


20 30 10 


SO AO 



9U 

• + 


311 

Oil 

on 

on 

□ ii 
on 
on 
oil 

on 

oil 

oil 

on 

on 

on 

on 

on 

on 

on 

□ n 

□ n 
on 


nnininnininn nnmnnnnnnnunnnnnnnnnnnnnnnnn 
n nnm in niinutininnninnnnnnnnnninntnnninnnnnn 
niiinum mini nnnnnnt i nun nun mini uni nnnnnnnnn n 
uuuuuuuiuiuuuuuiuuiuuninunuuuiu iinunnuuiuuunuiun 
luuuiuuinniuiuniunuuunuuunuunnuuiunu imiuummunii 
uuuunuumuinuuuuuuuuunuuuiuuiuuiunnnuuiuuuiiuuu 

iiuunuuuinuuuui 

uuuuuuiu juuuun 

nuuuuuuimuimii 

iimnn 

unnut 

uunm 


uiuuuuuniuuunuuunu 

uuuuunnuuuuuuuuun 

luuuinuuuuuuuuuuuu 


o 10 


20 30 10 


50 


AD 70 


aa 


90 


100 


1 

4 

♦ 

l 

* 

4 

i 

* 

4 

J 

4 

4 

I 

♦ 

4 

1 

4 

4 

J 

4 

4 

4*4 

I0u 


4-78 


SSTAT 

KM I N 

XH AX 

X0AR 

VAR 

SOEy 

iEHD 


STEP- ,147743*5 

CENTEttPOlNT OF INITIAL GROUP- ,07254016 

CENHRP01HT OF FINAL GROUP- 1.10*74572 


NO, OF 0B5£KVAUUnS« 3b 
NO, OF GROUPS- 7 


PERCENTAGE 

0 

10 

20 

30 


40 

50 

60 

70 

ao 

90 

too 

frequency 

0 

3 

7 

10 


14 

17 

21 

24 

28 

31 

35 

2*203-01 

l a »• * 



■••*Q{** 

• •••1 

It* 1 « **« 




••••• ( C»*«H 




4 

4 

4 

4 

X 

4* 

* 

4 

4 

4 

4 

4 


4 

4 

3* 

4 


♦ X 

4 - 

4 

4 

4 

4 

4 


4 

4 

♦ 

♦ 


+ 

x + 

♦ 

4 

4 

4 

4 


3*600-01 


5,158-01 


6.635-QI 


8*113-01 


9*540-0 1 


+ X 


4- 

x+ 

,*lo« 

* 


X* 

x+ 

1*0* 

♦ 

♦ 


+ 

>*!•* 

+ 

+ 

X* 

+x 

**l*l 

♦ 


♦ 

«$• 1 
♦ 
+ 
*■ 
♦ 


1. 107*00 | •••• 




1 •••* 






•t-aU-a 1 

FREQUENCY 0 

3 

7 

10 

14 

1 7 

21 

24 

2<S 

31 

35 

PERCENTAGE 0 

10 

20 

30 

40 

50 

60 

70 

80 

90 

100 


- .72S40159E-01 , 

- , I |O67457&*01 , 

- ,4t?M858E*DO, 

« .9H091H50E-0I , 

- ,3067433OE*tJO, 



a 15 THt ELEMENT ANALYZED 


sample 

INTENSITY 

ABSORPTION 

absorption 

AREA 

HEIGHT 

NUMBER 

COUNT 

FACTOR 

correction 

DENSITY 

PERCENTAGE 

1 

' 7 } 5 . 
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.99184 

• 9B2B4 

•94305-07 

• 17 

3 

1019. 
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